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FOREWORD 


Retardation  ayatems  uaed  £or  very  high  altitude  at  low  dynamic  pressurea  or 
for  low  altitude  cargo  delivery  often  require  exceptional  teat  equipment  and/or 
facilities.  As  tests  become  more  complex,  higher  developmental  costs  are 
incurred.  This  report  describes  a  technique  which  theoretically  permits 
parachute  systems  to  be  tested  under  more  convenient  conditions  of  altitude, 
velocity,  and  system  mass  resulting  in  the  same  conditions  of  opening  shock 
force  and  parachute  stress  distribution. 

Several  unique  properties  of  solid  cloth  parachutes  are  developed  along 
with  the  criteria  for  alternate  altitude  testing.  Methods  of  calculation  are 
provided  and  test  techniques  proposed. 
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INTRODUCTION 

\ 

The  flight  testing  of  solid  cloth  parachutes  sometimes  calls  for  unique 
test  conditions  such  as  very  high  altitudes  for  extremely  low  dynamic  pressures, 
or  heavy  payloads,  such  as  cargo,  tested  at  low  altitudes.  Each  of  these 
conditions  may  require  exceptional  test  equipment  and/or  facilities  such  as 
balloon-borne  and  rocket-assisted  high-altitude  vehicles  or  large  cargo  aircraft 
for  low-altitude  heavy  payloads. 

These  requirements  can  add  significantly  to  the  overall  test  program  cost 
and  complexity.  In  order  to  reduce  development  test  costs  and  complexity,  a 
theoretical  study  has  been  conducted  to  determine  the  possibility  of  altering 
the  test  conditions  in  a  manner  which  would  permit  the  use  of  lower  cost  test 
assets  and  still  obtain  the  same  p  irachute  opening  shock  forces  and  stress 
distributions.  The  results  of  this  study  show  that  if  the  system  mass  and 
test  velocity  are  varied,  as  shown  in  the  accompanying  method  for  a  constant 
Ballistic  Mass  Ratio  (BMK)  scale  parameter,  the  opening  shock  force  and  stress 
distributions  are  constant  for  all  altitudes. . . .  - 

The  examples  in  the  study  indicate  that  high-altitude  testing  may  be 
accomplished  at  low  altitudes  by  reducing  the  test  velocity  and  increasing  the 
payload  mass,  and  low-altitude  cargo  parachute  systems  may  be  tested  by  reducing 
the  payload  weight  and  increasing  the  test  velocity  and  altitude. 

In  the  case  of  the  modified  high  altitude  testing,  a  test  method  is 
proposed  which  improves  test  observation  and  control.  Modified  cargo  testing 
permits  several  test  vehicles  to  be  mounted  on  the  wing  racks,  of  an  A7  type 
aircraft  for  example,  with  properly  weighted  existing  test  vehicles  for  the  same 
total  weight  of  one  conventional  test  from  a  C-130  aircraft.  Eventually  full 
systems  tests  may  be  desired;  however,  parachute  development  time  and  costs  can 
be  reduced  by  alternate  altitude  testing. 

As  the  study  progressed,  several  unique  effects  of  altitude  on  solid  cloth 
parachute  performance  were  developed.  These  effects  are: 

a.  The  Ballistic  Mass  Ratio  as  a  scale  parameter. 

b.  Variation  of  opening  shock  force  at  variable  and 
constant  Ballistic  Mass  Ratios. 

c.  Parachute  stress  distributions. 
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d.  Average  steady-state  pressure  coefficients  and  the 
variation  as  a  function  of  the  inflation  time  ratio,  t/tQ, 

e.  Parachute  inflation  distance  and  inflation  reference  time. 

f.  System  impulse  as  a  function  of  Ballistic  Mass  Ratio. 

g.  Parachute  performance  as  related  to  the  inflation  time 
ratio,  t/to,  for  a  constant  Ballistic  Mass  Ratio. 

Although  the  study  was  specifically  for  solid  cloth  parachutes,  the  results 
may  be  considered  as  guideposts  for  trends  in  the  testing  of  other  types  of 
parachutes . 


APPROACH 


The  study  is  baaed  upon  two  conditions:  (a)  the  premise  that  the  Ballistic 
Mass  Ratio  satisfies  stated  criteria  required  for  a  genuine  scale  parameter,  and 
(b)  an  opening  shock  force  analysis  which  was  presented  at  the  fourth  AIAA 
Aerodynamic  Decelerator  and  Balloon  Technology  Symposium  in  1973.  Since  this 
paper  is  to  be  used  extensively,  it  is  included  as  Appendix  A  of  this  report  for 
the  convenience  of  the  reader.  Appendix  B  provides  a  worksheet  for  use  as  a 
guide  on  how  to  most  effectively  use  Appendix  A.  The  development  of  the  subject 
matter  is  to  be  a  combination  of  theoretical  development  combined  with  examples 
to  demonstrate  the  results.*  It  is  suggested  that  for  best  understanding  of 
this  report  the  reader  should  review  Appendix  A  first.  Essential  formulae  from 
Appendix  A  are  reviewed  below  as  a  basis  for  development. 


THE  BALLISTIC  MASS  RATIO  AS  A  SCALE  FACTOR 

A  method  of  scaling  parachute  inflation  and  steady-state  performance  has 
been  a  long-sought  element  of  parachute-system  analysis.  In  the  search  for  a 
scaling  parameter,  a  most  important  requirement  must  be  kept  in  mind.  That  is, 
that  for  any  quantity  to  be  a  valid  scaling  parameter,  it  must  relate  to  the 
variables  which  affect  performance. 

In  model  testing  of  airplanes,  missiles,  ships,  etc.,  in  wind  tunnels  and 
towing  basins,  Reynolds  number,  Froude  number,  and  Mach  number  are  acccepted 
scale  factors.  These  parameters,  together  with  ballistic  coefficient 
(W/CqSq)  and  surface  loading  (W/Sq),  have  not  provided  a  general 
correlation  of  parachute  test  data.  A  viable  scaling  parameter  for  general  use 
in  deployable  decelerator  testing  has  not  been  developed,  although  some 


*Where  formulae  from  Appendix  A  are  cited,  the  parenthesis  in  the  right  hand 
margin  of  the  page  denotes  the  appendix  page  and  formula  number,  respectively. 
Hence,  (A7-13)  reads  as  page  A7  formula  13.  Formulae  marked  with  a  single 
number  in  parenthesis  pertain  to  the  presented  development. 
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experlmencers  have  successfully  modeled  particular  cases.  A  look  at  the  methods 
of  testing  gives  us  a  clue  as  to  why  the  aforementioned  scale  factors  do  not 
apply  to  parachute  deployment.  In  the  testing  of  airplanes,  missiles,  ships, 
etc.,  a  rigid  model  is  constructed  and  mounted  in  the  test  facility  where  data 
is  recorded  at  one  or  more  constant  velocities.  Ttius,  three  parameters  which 
are  important  to  parachute  testing  are  eliminated.  First,  the  geometry  of  a 
deploying  parachute  undergoes  dramatic  changes  as  compared  to  the  constant 
geometry  of  a  rigid  model.  Second,  the  deploying  parachute  geometry  is  a 
function  of  deployment  time.  Third,  the  velocity  profile  obtained  during  a 
parachute  deployment  is  dependent  upon  the  mass  of  the  assembly  being  retarded. 
Transient  geometry,  time,  and  weight  are  not  factors  in  the  aforementioned  scale 
parameters.  This  leads  to  some  very  important  conclusions: 

a.  Reynolds  number,  Froude  number,  and  Mach  number  are  valid  scale  factors 
for  airplanes,  missiles,  ships,  etc.,  because  they  contain  the  variables  which 
affect  performance. 

b.  These  established  scale  parameters  are  unsuitable  for  parachute 
deployments  because  they  do  not  contain  important  variables  which  affect 
parachute  performance. 

c.  For  any  quantity  to  be  valid  scale  parameter  in  any  process,  it  must 
contain  those  variables  which  affect  system  performance. 

A  fundamental  example  of  a  “mass  ratio"  scale  factor  can  be  developed  by 
means  of  a  theoretical  horizontal  point  mass  trajectory.  For  simplicity,  only 
the  steady  state  (constant  drag  area,  C^Sq)  will  be  considered;  however,  the 
developed  mass  ratio  is  also  applicable  for  parachute  deployments. 

Example  \:  Determine  the  velocity  profile  of  an  automobile  which  is  being 
retarded  on  a  norizontal  road  by  a  fully  deployed  parachute.  Assume  negligible 
automobile  aerodynamic  drag  and  road  friction  forces. 

With  riference  to  Figure  1. 
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W 

FIGURE  1.  RETARDED  AUTOMOBILE 
OF  EXAMPLE  1 

The  quencicy  ZW/pgVgtCoSQ  is  in  effect  a  ratio  of  masses.  It  can  be 
visualized  as  in  Figure  2  to  be  the  ratio  of  the  mass  of  the  system,  W/g,  to  the 
mass  of  air  contained  in  a  right  circular  cylinder  of  face  area  Cj)So,  length 
Vgt,  and  mass  density  p.  If  this  quantity  is  denoted  by  M,  then  Equation 
(1)  becomes 


Vs 


(2) 


and  it  is  seen  that  the  velocity  profile  is  only  a  function  of  the  mass  ratio. 
Figure  3  illustrates  the  concept  of  how  the  air  mass  is  affected  as  the 
automobile  moves  along  the  retarded  trajectory.  As  time  increases,  additional 
air  mass  is  affected,  and  at  each  instant,  a  definite  velocity  ratio  and  mass 
ratio  exist,  as  shown  in  Equation  (2).  The  mass  ratio  contains  the  basic 
variables  (altitude,  p;  parachute  aerodynamic  size,  CpSp;  system  mass, 

W/g;  velocity,  V^;  and  time)  necessary  to  define  a  scale  factor  for  deployable 
decelerator  application.  Apparent  mass  is  another  effect  which  should  probably 
be  included,  but  is  not  presently  sufficiently  understood  or  defined  as  to 
permit  inclusion  in  the  analysis.  For  the  inflating  parachute  the  methods  for 
obtaining  the  dynamic  drag  area  ratio  signature  include  the  effects  of  apparent 
mass . 


Note  that  the  length  Vgt  in  Figure  3  is  not  a  true  trajectory  distance 
because  of  the  use  of  the  constant  initial  velocity  Vg .  The  true  distance 
(TD)  is: 
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FIGURE  2.  VISUALIZATION  OF  THE  MASS  RATIO  CONCEPT 


Xit  FIGURES.  VARIATIONOF  AFFECTED  AIR  MASS  ALONG  A 

FULLY  INFLATED  PARACHUTE  TRAJECTORY 
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Appendix  A  describes  a  mechud  of  calculating  the  opening  shock  forces  of 
several  types  of  solid  cloth  parachutes  which  utilizes  the  developed  mass  ratio 
as  applied  to  the  inflation  stage  of  operation.  Equations  (13)  and  (14)  on  page 
A-7  of  Appendix  A  relate  the  parachute  geometry,  air  flow  properties  of  the 
canopy  cloth,  drag  area  signature,  and  deployment  conditions  to  the  inflation 
reference  time  t^.  The  dependency  of  the  mass  ratio  on  the  system  parameters 
is  illustrated  in  the  flow  chart  of  Figure  4.  Therefore,  the  mass  ratio 
fulfills  Che  scaling  requirement  Chat  it  relates  to  Che  variables  which  affect 
performance  and  is  a  valid  scale  factor.  Since  the  mass  ratio  determines  the 
ballistic  performance  of  the  retarded  trajectory  the  name  "Ballistic  Mass  Ratio" 
(BMR)  is  appropriate. 

Although  Che  present  discussion  is  limited  to  solid  cloth  parachutes  the 
Ballistic  Mass  Ratio  approach  is  valid  for  all  types  of  parachutes.  Testing  at 
a  constant  BMR  at  all  altitudes  requires  a  modification  in  system  weight  to 
offset  changes  in  air  density  and  inflation  distance  (VgC^).  For  other 
types  of  parachutes,  Che  definition  of  Che  variation  of  inflation  distance  with 
altitude  is  Che  key  to  alternate  altitude  testing.  The  basic  idea  is  not 
complex.  Since  it  is  required  that  tests  are  to  be  conducted  at  a  constant 
dynamic  pressure,  Che  infIa;:ion  distance  associated  with  the  required  density 
and  velocity  is  defined.  This  causes  a  change  in  Che  BhJl.  The  system  mass  is 
then  adjusted  accordingly  to  return  the  BMR  to  the  original  value. 

When  a  given  solid  cloth  parachute  system  is  tested  at  constant  dynamic 
pressure  over  a  wide  range  of  increasing  altitudes,  a  decrease  in  inflation  time 
and  an  increase  in  opening  shock  force  occurs.  An  explanation  of  this  phenomena 
is  as  follows.  The  opening  shock  force  may  be  expressed  as 
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F  ^  F  X  * 

'max  's  *imax 

Fg  =  q  CdSq 

where  shall  be  shown  to  be  a  function  of  the  drag  area  ratio  at  line 

stretch,  n  =  CqS^/CjjSq,  and  the  BMR.  The  shock  factor  determines  the 
percentage  of  the  steady-state  drag  force  that  will  be  realized  as  opening  shock 
force.  For  altitudes  above  sea  level  there  is  a  decrease  in  air  density  and 
inflation  reference  time.  Both  of  these  effects  increase  the  BMR  and  the 
maximum  shock  force.  Testing  at  a  constant  BHR  results  in  the  same  opening 
shock  force  at  all  altitudes. 

The  decrease  in  inflation  reference  time  is  due  to  the  decrease  in  air 
density.  For  constant  dynamic  pressure  as  altitude  increases  the  test  velocity 
must  be  varied 


From  Equation  (14)  of  page  A-7  of  Appendix  A 


14W 

PgVgCoSo 


1 


or 


to 


14W 


eSVsL 


gpVo 

- 1 

n 

o 

c/> 

o 

_ 1 

2W 

/c  p\  1/2 

,Amo  ■  12“ j 

.e 

'  / 

-1. 

At  all  altitudes  the  Vg^  remains  constant,  and  the  primary  variable  is 
density.  Other  possible  variations  may  be  due  to  drag  coefficient  rise  caused 
by  a  decrease  in  canopy  cloth  air-flow  properties,  and  possible  changes  in 
geometry  caused  by  the  rise  in  drag  coefficient.  Otherwise,  the  reference  time 
is  only  a  function  of  system  constants  and  altitude. 


DISCUSSION  OF  APPLICABLE  APPENDIX  A  FORMULAE 

During  the  inflation  of  a  parachute  canopy  the  instantaneous  force  can  be 
described  as  the  steady  state  drag  force,  Fg ,  times  some  instantaneous  shock 
factor,  x^j^. 


F  -  FgXi 
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where  X£  at  any  instant  is  a  function  of  the  instantaneous  drag  area 
signature,  CqS/CqSq,  and  the  instantaneous  velocity  ratio  (V/Vg)^. 

The  drag-area  signature  is  indicative  of  the  type  of  parachute  being  used  and  is 
dependent  on  the  deployment  time  ratio,  t/tg,  and  the  method  of  deploying  the 
parachute  described  by  the  initial  drag-area  ratio,  n  "  CqS^/CqS^. 

For  the  types  of  solid  cloth  parachutes  of  Appendix  A  the  general  case  drag-area 
signature  was  determined  to  be 


CpS 

CpSo 


2 

(AA-4) 


The  parachute  drag-area  signature  may  be  presently  obtained  by  two  methods: 

a)  Infinite  mass  wind-tunnel  tests  wnere  the  methods  described  on  pages 
A-3  and  A-4  are  used. 

b)  Free  flight  tests  where  event  times  have  been  established  and 
simultaneous  measurements  of  parachute  force  and  dynamic  pressure  have 
been  recorded  during  deployment. 

Tnen: 


CpSx 


F<p 

‘‘t 


and  Che  drag  area  signature  at  any  instant 


CpSx  Fx  CpS 

CpSo  *lx^D^o  ^D^o 


This  writer  prefers  the  wind  tunnel  rethod  as  it  offers  better  control  of  test 
conditions  and  events,  but  the  second  method  may  offer  a  better  utilization  of 
existing  field  test  data  over  a  range  of  altitudes,  velocities,  and  various 
types  of  parachutes. 

The  importance  of  Che  discovery  that  the  dynamic  drag-area  signature  and 
geometry,  as  a  function  of  the  time  ratio,  c/Cq,  are  independent  of  altitude, 
velocity,  and  system  mass  has  not  yet  been  fully  realized.  This  independence 
separates  Che  geometry  of  deployment  from  the  forces  generated  during  inflation, 
and  also  Che  stress  distributions  in  the  canopy. 
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For  simplicity  in  demonstrating  the  principles  of  this  analysis  n  is  to  be 
considered  as  zero. 

Hence , 


Cd^o  V^o/ 

Utilizing  Newton's  third  law  of  motion  in  the  horizontal  test  attitude*  results 
in: 


L-  ^M-dt _ ^ _ V  ^ 

to  ^  PgVgtoCDSo  •  v2 

which  leads  to  the  velocity  ratio  equation: 


V_  _  _ 1 


(A2-2) 


(A4-7) 


and  the  shock  factor: 


*i 


(A5-8) 


Equation  (A2-2)  also  utilized  what  has  come  to  be  known  as  a  Ballistic  Mass 
Ratio  (BMR). 


M 


2W 

PgCoSoV.to 


(A4-6) 


*Parachutes  tested  in  vertical  fall  exhibit  an  opening  shock  force  of  Ig 
greater  than  horizontally  deployed  canopies.  This  is  a  limitation  on  the 
development . 
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One  of  the  verieblee  in  the  BMR  is  the  inflation  distance,  VgtQ,  which  is 
the  product  of  the  system  velocity  at  line  stretch,  V^,  and  the  reference 
opening  time,  t^. 


8PYo 

CrSo 

2U 

Amo  ■  Aso  2  / 

V  t 

*  ®  PgCoSo 


(A7-14) 


Note  that  when: 


0 


^Ho  “  ^o  ^ 


1/2 


(4) 


a  critical  inflation  condition  exists,  and  at  any  test  velocity  the  t^  time  is 
infinite.  As  altitude  increases,  the  density  p  approaches  zero  and  equation 
(4)  becomes  positive  and  the  parachute  will  inflate.  For  imporous  canopies  k^O 
and  the  parachute  always  inflates.  These  boundary  conditions  agree  with  general 
field  test  experience. 

The  inflation  distance  is  not  dependent  on  the  test  velocity  Vg.  Rather  it  is 
a  function  of  system  weight,  w,  and  density  altitude,  0;  the  steady-state 
geometric  characteristics  of  the  particular  type  of  parachute,  and  drag  area, 
CpS^;  voluate,  nouth  area,  A{(o;  Surface  area,  Ago'Sp;  cloth 
airflow  properties,  k  and  n;  and  average  pressure  coefficient,  Cp,y< 

The  average  pressure  coefficient  for  the  steady-state  parachute  was  determined 
as  follows: 

The  steady-state  drag  force,  in  a  wind  tunnel  for  example,  may  be  written 

but  it  may  also  be  written 
F  -  AP«vSp 

where  Sp  is  the  projected  area  of  the  inflated  canopy  and  AP  is  the  average 
pressure  differencial  acting  on  8p* 


AP,vSp  ■  F  ■  qCpSo 


AP.v  CpSo 
^I'av  q  ”  Sp 
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Ratios  of  Dp/Dg  from  Reference  2  are  presented  in  Table  1  for  several  types 
of  parachutes. 


TABLE  1.  RANGE  OF  AVERAGE  STEADY  STATE  CANOPY  PRESSURE  COEFFICIENTS 


PARACHUTE 

TYPE 

Co 

RANGE 

t>p/Oo 

RANGE 

Cpt» 

RANGE 

FLAT 

0.75 

0.70 

1.631 

CIRCULAR 

0.80 

0.87 

1.782 

EXTENDED 

0.78 

0.70 

1.692 

SKIRT 

0.87 

0.66 

1.997 

CROSS 

0.60 

0.72 

1.167 

0.78 

0.66 

1.791 
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For  a  flat  circular  parachute  with  a  C[3f*0.75  and  Dp/Do"0.67 

r  -  0-75 

(0.67)2 


Cp  -  1.671 

It  may  be  concluded  that  the  average  steady-state  canopy  pressure 
coefficient  is  the  canopy  drag  coefficient  based  on  the  projected  area.  The 
inflation  geometry  as  a  function  of  t/t^  is  constant  with  altitude.  It  can  be 
deduced  from  this  that  if  Che  inflating  geometry  is  independent  of  altitude,  Che 
forces  which  cause  the  geometry  Co  develop  must  also  be  independent  of 
altitude.  These  forces  come  from  the  pressure  distribution  along  Che  gore 
panel.  Therefore,  at  any  given  c/c^  a  definite,  repeatable  pressure 
distribution  exists  which  varies  with  t/t^,  but  is  constant  for  all  altitudes 
and  Che  is  constant  with  altitude. 

The  Air  Force  Plight  Dynamics  Laboratory  conducted  field  tests  of  a  28-fooc 
(Dq)  Solid  Flat  Circular  Parachute  (SFCP)  system  at  altitudes  of  6,000,  13,000, 
and  21,000  feet.  A  result  of  these  tests  was  the  demonstration  of  the 
repeatability  of  Che  canopy  projected  area  to  surface  area  ratio  at  all  test 
altitudes.  With  the  data  of  Figure  5  the  average  steady-state  pressure 
coefficients  for  the  28-fooc  SFCP  can  be  estimated. 
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FIGURE  6.  NORMALIZED  CANOPY  AREA  GROWTH  DURING  INFLATION  OF  28  FT  ID^)  SOLID  FLAT 
CIRCULAR  PARACHUTE 
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Cp 


av 


\3 

2 

Cd 

a-„)  (f- 

/  ^  " 
3' 

t/tT 

0.0117  +  0.0034  (118.4)  ° 


(6b) 


Figure  6  shows  the  estinated  Cpi^y  for  the  following  conditions:  Cq  ■  0.75; 
n  ■  0.0024;  tg  ■  t£;  and  the  two  variations  of  Sp/S^  as  a  function  of 
t/t£  in  Figure  3.  Cloth  air  flow  properties  "k"  and  "n"  may  be  evaluated  by 
the  methods  of  formulae  (A13-29)  and  (A13-30).  Another  method  is  to  adapt  a 
least  squares  fit  through  the  data  points.  A  true  average  value  of  k  sad  n  in 
any  given  parachute  requires  a  large  number  of  permeameter  rate  of  air  flow  tests 
due  to  the  variation  of  the  woven  cloth  air  flow  rates.  Their  values  may  also 
vary  for  different  parachutes  in  the  same  lot.  This  is  most  likely  one  of  the 
causes  of  performance  variations  for  similar  test  conditions. 


FIGURE  6.  ESTIMATED  STEADY  STATE  AVERAGE  PRESSURE  COEFFICIENT  FOR  THE 
INFLATING  28  FT  (Do)  SOLID  FLAT  CIRCULAR  PARACHUTE  OF  FIGURE  5 
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The  inflation  time  and  distance  may  also  be  calculated  by  equation 
(A7-13).  In  this  equation  the  actual  value  of  "n"  for  cloth  airflow  prooerties 
is  used  rather  than  the  assumed  value  of  '*n''«l/2.  At  low  altitudes  (A7-13) 
yields  a  better  value  of  inflation  distance  than  (A7-14). 


DEVELOPMENT  OF  ALTERNATE  ALTITUDE  TESTING 

The  drag-area  signature  is  independent  of  altitude,  velocity,  and  system 
mass,  Reference  1,  while  the  velocity  ratio  is  dependent  on  the  BHR,M.  Since 
the  shock  factor  is  a  function  of  the  velocity  ratio,  it  is  also  a  function  of 
the  BMR,M.  For  alternative  altitude  testing,  the  forces  generated  may  be 
duplicated  at  other  density  altitudes  by  testing  at  a  constant  BMR.  The 
requirements  for  generating  a  constant  BMR  are  to  adjust  the  system  weight  for 
the  particular  density  altitude  desired.  Note  that  density  altitude  is  stressed 
as  this  may  vary  from  geometric  altitude  or  pressure  altitude  used  for  some 
altimeters . 

Replacement  of  the  inflation  distance  in  the  BMR  equation  (A4-6}  by  its 
equivalent  (A7-14)  results  in: 


BPVo 

CdSo 

■ 

2W 

_Amo  - 

> 

CD 

o 

y 

p)'" 

_ 2W _ 1 

14W  K1  ,1  7  TrKl-n 

‘  L' -‘J  L  !l 
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K1 


1 

1  7M 


K1  “ 


(7) 


Solving  for  the  required  weight  as  a  function  of  density  altitude. 


(8) 


The  inflation  distance  at  any  alternate  density  altitude  is  determined  by 
equation  (A7-14)  and  the  required  weight. 


VsCo 


14W 

PgCoSo 
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(9) 


DETERMINATION  OF  REQUIRED  TEST  VELOCITY 

Since  the  steady-state  drag  force  is  a  function  of  altitude,  the  steady- 
state  force  will  remain  constant  if  the  tests  are  conducted  at  constant  dynamic 
pressure.  The  dynamic  pressure  at  the  original  density  altitude  must  be  equal 
to  the  test  altitude. 

original  qg  *  test  qgj 


VsT  "  Vg 

p/ PSL  *  density  ratio  at  original  requirement  altitude 
PT/psl  *  density  ratio  at  alternate  altitude 


(10) 


INFLATION  REFERENCE  TIME 

The  inflation  reference  time,  t^,  at  each  altitude  is  calculated  for  the 
conditions  at  that  altitude. 


Example  2:  In  order  to  illustrate  the  theory  developed,  a  sample  problem 
shall  be  used. 

Problem:  A  T-10  type  parachute  is  to  be  tested  horizontally  at  a  density 
altitude  of  80,000  feet  at  300  feet  per  second  test  velocity.  The  system  weight 
is  250  Ib. 

Determine:  Alternate  test  altitudes  to  obtain  the  same  maximum  opening  shock 
force . 

Solution:  For  a  T-IO  type  parachute  the  diameter  is  D„*35  ft  and  the  drag 
area  is  721  ft2 

Determine  the  steady-state  parachute  geometry,  Opg^ ,  and  cloth  air-flow 
constants  k  and  n. 

From  Table  2,  page  (A14)  for  a  10  percent  Extended  Skirt  parachute  of  30  gores 


17 


NSWC  TR  85-2A 


-  0.650  ;  J  *  0.825  ;  |  -  0.6255  ;  ~  -  0.7962 


S  ,  ,  0^65^  35  .  ^,^3,3 


The  steady-state  mouth  area, 


Amo  •  t 


/n/I  -  b/aV 

^  ■  \  b'/i  / 


Afio  -  i(ll.375)‘ 
Amo  -  381  ft2 


(0.825  -  O.6255Y 
0.7962  \ 


(12) 


The  steady-state  canopy  volume, 


Vo  -  2  ^ 


(A14-31) 


Vq  -  I  »(11. 375)3  (0.6255  +  0.7962] 


-  4382  ft3 


The  cloth  air-flow  constants,  k  and  n,  were  determined  for  the  Mil-C-7020,  type 
III  Cloth,  using  the  methods  on  page  A-13. 


87.6255  ,  ,, 

k  - 1.46 


n  -  0.63246 


The  term,  k,  as  determined  for  the  nominal  permeability  has  units  of  ft/min. 
However,  calculations  are  based  on  ft/sec.,  hence  the  division  by  60.  The 
average  pressure  coefficient  was  taken  as 

Determine  the  BMR  at  80,000  feet 

a.  Use  equation  (A7-14)  to  determine  inflation  distance 
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The  results  of  these  calculations  are  summarized  in  Table  2  and  graphically 
illustrated  in  Figure  7. 

The  solution  of  example  2  using  eq.  (A7-13)  requires  the  use  of  a  computer 
due  to  the  value  of  "n". 


Vo 


%„VstoM  ln[l  ♦  -  A,o  f 


dt 


[‘  *  I.  (i)T 


A  method  of  calculation  is  as  follows: 


(A7-13) 


a)  For  the  given  density  altitude  determine  the  test  velocity,  Vgj, 
eq  (10). 

b)  Let  dt  =  to/10,000 

c)  Assume  a  value  of  t^  and  compute  ^  for  the  system  parameters 
listed  in  the  heading  of  Table  2. 

d)  Compare  the  computed  canopy  volume  with  the  previously  calculated 
canopy  volume  of  4382  ft^. 

e)  If  the  computed  canopy  volume  does  not  agree  with  the  4382  ft^  within 
a  specified  limit  of  5  ft^  then  correct  to  end  iterate  again. 

f)  Let 


4382  t 

o 

to  »  - - 

•^o  computed 

g)  Calculate  the  inflation  distance. 

Vgto  -  Vg  X  to 

h)  Determine  the  weight  required  for  M  =  2.950 


W 


pgC_S  V  t  M 
°  D  o  s  o 

2 


The  integral  of  equation  (A7-13)  has  an  interesting  property  in  that  if  the 
integrated  function  for  any  altitude  is  divided  by  the  to  time  for  that 
altitude  a  constant  value  at  all  altitudes  is  obtained,  as  shown  in  Table  3. 
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WEIGHT  -  250  LB 
Dq  =  35  FT 

-  S„  -  962  FT2 


DENSITY 

ALTITUDE 

KILOFEET 


0.03606 


0.05856 


0.09492 


.15311 


0.24708 


0.37473 


0.53316 


0.73859 


00 


TEST  CONDITIONS  @  60,000  FEET 

A^q"381FT^  Ic-1.46 

C-^  =  721  FT^  n  0.63246 

T2^-4382Ft3  Cp,^'1-7 


Eq  (A6-  14);n  »  1/2 


VsT 

W 

FPS 

LB 

300.0 

250.0 

235.4 

409.7 

184.9 

671.6 

145.6 

1099.4 

114.6 

1808.0 

93.1 

2800.3 

78.0 

4072.9 

66.3 

5783.1 

57.0 

8048.6 

M  -  2.950 
Vj  -  300  FPS 
Fg  -  2782  LB 


Eq  (A6-  13);  n  '•0.63246 


TABLE  3.  RATIO  OF  THE  INTEGRAL  OF  EQUATION  (A6-13)  TO  THE 
INFLATION  REFERENCE  TIME  FOR  EXAMPLE  2 
✓  t- 


r  - ^ 

’"VvJ 


M  >  ^950 
n  -  0.63246 


ALTITUDE 
THOUSANDS 
OF  FEET 


0.28713 

0.37048 

0.47924 

0.62114 

0.81041 

1.02758 

1.26578 

1.54312 

1.86861 


0.04256 

0.05492 

0.07104 

0.09207 

0.12013 

0.15232 

0.18762 

0.22874 

0.27608 


0.14823 

0.14823 

0.14823 

0.14823 

0.14823 

0.14823 

0.14823 

0.14823 

0.14823 


0.63246 


FIGURE  7.  EFFECT  OF  CLOTH  AIR  FLOW  PARAMETER  "n”  ON  THE  REQUIRED  SYSTEM  WEIGHT.  INFLATION  DISTANCE.  AND  INFLATION 
REFERENCE  TIME  FOR  ALTERNATE  DENSITY  ALTITUDES  AT  CONSTANT  BALLISTIC  MASS  RATIO  FOR  EXAMPLE  2 
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The  integral  of  "B"  is  illustrated  in  Table  4  for  various  values  of  BHR  and  n« 
With  these  data  the  t^  time  of  eq  (A7-13)  may  be  calculated. 


^o  “  *^o 


Atio  Vs  M  In  [l  >  y  -  Ago  ^8^“  (|^) 


For  an  imporous  solid  cloth  parachute  canopy,  k  °  0  and 


(13) 


^o 


V 

-o 


,  r 

- 1' 

^Mo  Vg  M  In 

(13a) 


Assume  a  value  of  t^  and  calculate  the  BMR.  For  the  known  values  of  “BHR”  and 
"n"  interpolate  "B/t^"  from  Table  4  and  calculate  t^  using  eq  (13).  Figure  8 
illustrates  a  method  for  limiting  the  number  of  required  calculations  by 
plotting  the  assumed  and  calculated  t^  values.  The  desired  tg  is  the 
particular  value  at  the  intersection  of  Che  plotted  data  and  the  locus  of  Che 
Cq  assumed  equals  Cq  calculated  line. 


OPENING  SHOCK  FORCE 

If  Che  maximum  shock  force  occurs  in  the  finite  mass  regime,  then  the  time 
of  occurrence  during  the  inflation  process  is  determined  by; 


1/7 


and  Che  maximum  shock  factor 


’‘iraax 


(A3- 10) 


The  limiting  BMR  for  finite  mass  operation,  see  Appendix  C,  is  the  mass  ratio 
which  causes  the  maximum  shock  to  occur  when  t^tj,. 
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TABLE  4.  RATIO  OF  THE  INTEGRAL  OF  EQUATION  (A7-13)  TO  THE  INFLATION 
REFERENCE  TIME  FOR  VARIOUS  VALUES  OF  “BMR"  AND  "n" 


fl 

OAOOO 

OAOdO 

0.7000 

UJKOO 

O.SOOO 

1.0000 

0.  010 

0.  03737 

0.  02102 

0.  0l**0 

0.  01343 

0.  01106 

0.  00624 

0.  020 

0.  04165 

0.  03434 

0.  03039 

0.  02364 

0.  02023 

0.  01754 

0.  030 

0.  05355 

0.  04433 

0.  03770 

0.  03393 

0.  0303* 

0.  02476 

0.  040 

0.  04061 

0.  05243 

0.  04554 

0.  03669 

0.  03910 

0.  03129 

0.  050 

0,  04751 

0.  05614 

0.  05214 

0.  04427 

0.  04126 

0.  03704 

0.  040 

0.  07311 

0.  04468 

0.  05764 

0.  05164 

0.  04670 

0.  04226 

0.070 

0.  077B7 

0.  04662 

0.  04265 

0.  05461 

0.  05156 

0.  04660 

0.  060 

0.  06166 

0.  07413 

0.  04724 

0.  0*134 

0.  05595 

0.  05126 

0,  0^0 

0.  06556 

0.  07763 

0.  07116 

0.  04522 

0.  05962 

0.  05522 

0.  too 

0.  06674 

0.  06133 

0.  07471 

0.  0*001 

0.  06353 

0.  05662 

0.  125 

0.  09930 

0.  06636 

0.  00314 

0.  07446 

0.  07134 

0.  04*60 

0.  150 

0.  10036 

0.  06364 

0.  00000 

0.  003*0 

0.  07773 

0.  07216 

0.  175 

Oi  10446 

0.  06649 

0.  C9264 

0.  06761 

0.  06207 

0.  07067 

0.  200 

0.  10764 

0.  10223 

0.  06700 

0.  06213 

0.  06739 

0,  06236 

0.  225 

0  11045 

0.  10539 

0.  10049 

0.  06562 

0.  09146 

0.  06740 

0.  250 

0.  11304 

0.  10606 

0,  10341 

0.  06600 

0.  09465 

0.  09094 

0.  379 

0.  11510 

0  1 1041 

O.  10566 

0.  10176 

0.  06761 

0.  09405 

0.  300 

0.  11466 

0.  1 1245 

0.  10034 

0.  10423 

0.  10042 

0.  09661 

0.  339 

0.  1 1644 

0.  11425 

0.  11023 

0.  10*40 

0.  10275 

0.  09937 

0.  350 

0.  11665 

0.  11904 

0.  11301 

0.  10034 

0.  10404 

0.  10149 

0.  375 

0.  12110 

0.  1  1727 

0.  11340 

0.  11000 

0.  10672 

0,  10349 

0.  400 

0.  12221 

0.  11655 

0.  II9A4 

0.  11144 

0.  10642 

0.  10520 

0.  425 

0.  12321 

0.  1  1671 

0.  11434 

0.  11309 

0.  10667 

0. 10696 

0.  450 

0.  12413 

0.  12074 

0.  11793 

0.  11440 

0.  11130 

0.  10648 

0.  475 

0.  12464 

0.  12173 

0.  11641 

0.  11999 

0.  113*0 

0. 1C907 

0.  50C 

0.  13973 

0.  12241 

0.  1 19*0 

0.  114*9 

0.  11306 

0.  11116 

0.  939 

0.  13443 

0.  13343 

0.  13093 

0,  11771 

0.  11496 

0.  11235 

0.  990 

0.  13704 

0.  12417 

0,  12134 

0  110*9 

0.  11*01 

0. 11245 

0.  575 

0.  13744 

0.  13404 

0.  13319 

0,  11993 

0.  1 1*9* 

0.  11**0 

0.  400 

0.  I3«3I 

0,  13991 

0.  13300 

0.  13033 

0.  11709 

0.  119*3 

0.  625 

0.  12672 

0.  13411 

0,  12354 

0.  13100 

0.  110*7 

0.  11*33 

0.  450 

0,  13930 

0.  12447 

0.  13430 

0.  13179 

0.  11V49 

0, 11717 

0.  475 

0.  12645 

0.  12719 

0.  12476 

0,  13349 

0.  12016 

0. (1796 

0.  700 

0.  13007 

0.  12746 

0.  13939 

0.  12306 

0.  1300* 

0. 11670 

0.  739 

0.  13044 

0.  12614 

0.  12566 

0.  133** 

0. 12150 

0.  11929 

0.  750 

0.  I1003 

0.  13097 

0.  13*37 

0.  13433 

0.  13311 

0.  12005 

0.  775 

0.  13110 

0.  13090 

0.  12464 

0.  13474 

0.  133*0 

0. 12066 

0«  600 

0.  13.91 

0.  13937 

0.  13730 

0.  13933 

0.  13333 

0.  13137 

0.  625 

0.  13162 

0.  12674 

0.  12770 

0.  13970 

0.  13374 

0.  13193 

C.  050 

0.  13212 

0.  13006 

0.  12606 

0.  12414 

0.  12422 

0. 12236 

0.  675 

0.  13236 

0.  13041 

0.  12647 

0.  12454 

0.  12466 

0.  12266 

0.  900 

0.  13244 

0.  13072 

0.  12662 

0.  12464 

0.  18514 

0.  12234 

0.  423 

0.  13261 

0.  13103 

0.  12614 

0.  12724 

0,  12556 

0.  12360 

0.  950 

0.  13315 

0.  13130 

0.  12646 

0.  12771 

0.  12566 

0.  12424 

0.  975 

0.  13336 

0.  13197 

0.  12960 

0.  12605 

0.  12634 

0.  12466 

000 

0.  13340 

0.  13103 

0.  13006 

0.  13039 

0.  12671 

0.  12306 

i.  500 

0.  13453 

0.  13529 

0.  13407 

0.  13307 

0.  i3i*a 

0,  13090 

3.  000 

0.  13604 

0.  13711 

0.  13417 

0,  13934 

0.  13438 

0.  13240 

2.  500 

0.  13600 

0.  13623 

0.  13747 

0.  12471 

0.  13999 

0.  12521 

3.  000 

0.  13644 

0.  13699 

0.  13634 

0,  12770 

0, 13707 

0.  12644 

3.  500 

0,  14010 

0.  13954 

0.  13666 

0.  1264.1 

0. 13760 

0.  12723 

4.  OOO 

0.  14044 

0.  13695 

0.  13644 

0.  12667 

0.  12649 

0.  13600 

4.  500 

0.  14071 

C.  14027 

0.  13664 

0.  12640 

0.  12067 

0.  12654 

5.  000 

0.  14093 

0.  14054 

0.  14014 

0.  12675 

0.  13625 

0.  13696 

5,  500 

0.  14111 

0.  14075 

0.  14039 

0.  14003 

0.  13667 

0.  12928 

6.  000 

0.  14124 

0.  14063 

0.  140*0 

0.  14037 

0.  13994 

0.  12961 

6.  500 

0. 14139 

n,  14106 

0.  14077 

0.  14047 

0.  1401* 

0,  1399* 

7,  000 

0. 14150 

0,  14131 

0.  14063 

0.  14044 

0.  1403* 

0.  14006 

7.  500 

0. 14156 

0.  14133 

0.  14104 

0.  14076 

0.  14052 

0.  14026 

6.  000 

0.  14146 

0.  14143 

0.  14117 

0.  14092 

C.  14066 

0.  14042 

6.  500 

0.  14175 

0.  14151 

0.  14120 

0. 14104 

0.  14061 

0. 14057 

9.  000 

0.  14161 

0.  14156 

0.  14137 

0.  14115 

0.  14062 

0.  14070 

9  500 

0.  14167 

0.  14144 

0.  14145 

0.  14124 

0. 14102 

0. 14062 

10.  000 

0.  14193 

0.  14172 

0.  14152 

0.  14133 

0. 14112 

0.  14092 

10. 500 

0.  M1V7 

0.  14170 

0.  14156 

0.  14140 

0.  I4I3I 

0.  14102 

1 1.  000 

0,  14202 

0,  14163 

0.  S4I45 

0.  14147 

0.  14126 

0.  14110 

1 1  500 

0,  14204 

0,  14106 

0.  14171 

0.  14152 

0.  14124 

0.  14116 

12. 000 

0.  14206 

0.  14162 

0.  14174 

0,  14159 

0.  14142 

0.  7  4125 

12.  500 

0.  14213 

0,  14194 

0.  14t60 

0,  14144 

0. 14146 

0. 14132 

1 3.  000 

0.  14214 

0.  14200 

0.  14165 

0.  14166 

0.  14134 

0.  14136 

13.  500 

0.  14216 

0.  14204 

0.  14166 

0. 14174 

0. 14156 

0. 14144 

14.  000 

0.  14221 

0.  14207 

0.  14162 

0.  14176 

0.  14164 

0. 14149 

14.  500 

0.  14224 

0,  14210 

0.  14164 

0.  14162 

0. 14160 

0. 14154 

IS.  000 

0.  14226 

0,  14213 

0.  14166 

0.  14104 

0.  14172 

0.  14159 

LISTED  VALUES  ARE  VALID  FOR  ALL  ALTITUDES 
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ASSUMED  (SEC) 


FIGURE  8.  METHOD  FOR  LIMITING  THE  ITERATIONS  IN  CALCULATING 
BY  EQUATION  13 


Ml  •  ^  ■  0.1905  (for  n  “  o) 

BHR's  greater  Chan  Ml  require  chat  Che  canopy's  elascicicy  and  structural 
strength  be  taken  into  account.  The  opening  shock  is  defined  again  as 
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(A9-19) 


Where  the  BMR  is  the  same  as  before  and  the  subscript  “o”  refers  to  conditions 
at  the  time  t^.  The  point  here  is  that  the  shock  factor  is  still  a  function 
of  the  BMR  and,  that  for  a  constant  BMR,  the  shock  factors  and  velocity  ratios 
are  the  same  at  all  density  altitudes  throughout  the  inflation  process. 


The  maximum  shock  factor  in  the  elastic  range  o f  operation  is  determined  by 
the  canopy  constructed  strength,  and  the  ultimate  material  elongation, 

^max'  ^an  these  properties  remain  constant,  the  maximum  opening  shock  is 
constant  with  altitude. 


^imax 


fe)‘ 


.  4:  m  -I 


and  the  maximum  shock  force  is 


*^max  *  ^8  Hmax 

The  BMR  of  2.950  exceeds  the  finite  mass  limiting  mass  ratio,  Mi^«0.l905. 
Therefore  the  maximum  opening  shock  force  will  occur  after  time  t^,  and  the 
constructed  radial  strength  of  the  canopy  and  the  ultiisate  material  elongation 
must  be  included. 


Calculate  velocity  ratio  and  shock  factor  at  t'C^. 


Vo  I 

'  I  .  ^ 
7M 


0.9538 


(A9-18) 


(a’  ■  - 


(AlO-22) 


The  steady-state  drag  force  of  the  system  at  80,000  feet 


Fs  -  7  pVs^CdSj 
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Pg  -  ~  X  8.575  X  10"5  X  3002  x  721 
Pg  -  2782  lb. 


laitial  elongacioa  of  radials 


Mi 

*o  “  f  *max 


0.9098  X  2782 
“  30  X  550  * 


€o  »  0.0337 


Prom  Pigure  15  on  page  (AlO) 


where  P^ 

0.22 


30 


gores  X 


550  lb 
gore 


^pSmax 

CdSo 


1.07 


The  racio  of  inflacion  time,  Cf,  Co  Cg  is 

if  /cdS...Y^‘ 

to  ^  CdSo  ) 

tf  , ,, 

7”  -  (1.07)1/^ 

*•0 

—  -  1.01134 

^o 

The  maxiffluffl  shock  factor  aC  tf 


(AlO-24) 


(A9-19) 
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x; 


1.07 


Linax 


:^T3F"  *  TTI 


^.956 


- 7/6  -112 

07)  -1 


’^imax  "  0.9661 


The  aiaximum  shock  force  is 

^max  “  ^s  *imax 
^aax  ■  2782  x  0.9661 
^max  “  2688  lb 

COMPARISON  OF  OPENING  SHOCK  FORCES  AT  80,000  FEET  AND  AT  SEA  LEVEL 

It  has  been  shown  Chat  the  velocity  ratios  and  shock  factors  are  a  function 
of  the  BMR,  canopy  constructed  strength,  and  ultimate  material  elongation.  If 
these  quantities  remain  constant  for  all  alternate  altitudes  and  the  tests  are 
conducted  at  constant  dynamic  pressure,  the  sea  level  opening  shock  force  is 
exactly  equal  to  the  shock  force  at  80,000  feet.  The  inflation  times  over  which 
Che  forces  are  generated  vary  widely,  but  the  maxima  are  the  same.  This  effect 
is  illustrated  in  Figure  9.  When  the  drag  area,  velocity  and  dynamic  pressure 
ratios,  shock  factor,  and  force  are  considered  as  a  function  of  the  time  ratio, 
c/t0,  the  performance  during  inflation  is  identical  for  all  altitudes  as  shown 
in  Figure  10. 


EFFECT  OF  TEST  METHOD  ON  CANOPY  STRESS  DISTRIBUTION 

A  truly  effective  alternate  altitude  test  technique  should  provide  a  stress 
distribution  during  inflation  which  is  Che  same  as  Che  stresses  at  the  original 
altitude.  The  drag-area  signature  is  associated  with  a  definite  inflated 
geometry  at  each  t/tg  during  inflation.  Since  the  drag-area  signature  is 
independent  of  altitude,  the  associated  geometry  and  pressure  distribution  are 
likewise  independent.  Use  of  a  constant  BMR  at  alternate  altitudes  means  that 
Che  instantaneous  dynamic  pressure  ratio  is  the  same  at  all  altitudes.  At  a 
given  c/t0,  CdS/CdSq  is  constant  and  therefore  Che  geometry,  r/r^,  is 
constant,  and  q/q^  is  constant. 

o  ■  APr 


o  “  Cp  qr 
av 


FORCE  LBS 
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The  is  coasCanC  with  altitude,  then  at  a  given  t/tQ, 


(t) 

ik) 


f§2_.  1 

oalt 

The  stress  distributions  are  identical  for  all  altitudes. 


ALTERNATIVE  ALTITUDE  TEST  METHC^ 

Reduction  of  the  required  test  altitude  is  accomplished  by  an  increase  in 
the  system  mass  and  lowering  the  test  velocity  to  maintain  constant  dynamic 
pressure.  In  example  2  the  test  goals  may  be  achieved  near  sea  level  conditions 
if  a  test  weight  of  8048.6  pounds  and  a  test  velocity  of  57  fps  (38.9  mph)  were 
used.  The  problem  is  to  find  a  test  vehicle  which  can  lift  that  much  weight  at 
such  a  low  velocity.  Helicopters  are  one  possibility;  however,  Che  effects  of 
the  rotors  may  affect  test  performance.  Another  possibility  is  the  use  of  a 
truck  with  a  tower  mounted  on  it  ^rom  which  Che  test  mass  could  be  pulled  during 
parachute  deployment.  For  a  parachute  test  altitude  of  one  from  the  lower 
hem  of  the  parachute  Co  Che  ground,  a  cower  with  a  parachute  centerline  height  of 

Dq  -«■  a  s  46.4  feet  above  Che  ground  is  required.  A  test  rig  like  this  may  not 
be  the  most  desirable  system  from  Che  standpoint  of  personnel  safety  and  tower 
wake  effects.  A  third  possibility  is  the  use  of  a  sled  track  where  the  test 
assembly  can  be  launched  from  the  edge  of  a  cliff  into  undisturbed  air. 

Barometric  and  temperature  sensors  permanently  placed  in  the  test  area  can 
provide  a  continuous  monitoring  of  the  test  density  altitude  for  correct  vehicle 
weight  and  test  velocity.  The  test  vehicle  can  be  built  slightly  under  the 
programmed  test  weight  and  ballasted  prior  to  the  test  for  the  proper  mass  for 
the  teat  density  altitude.  The  test  mass  is  accelerated  to  a  planned  velocity 
by  a  truck  type  o f  reusable  pusher.  Upon  braking  the  pusher  disengages  from  the 
test  mass  which  slides  off  the  track  muzzle  at  test  velocity.  A  possible 
scenario  is  shown  in  Figure  11. 

At  sea  level  the  equilibrium  velocity  of  Che  system  is: 


Ve  - 

Ve  - 


jrz 

T  p  Cj)S 


002378 


8048.6 

721 
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FIGURE  It.  OFF  THE  CLIFF  METHOD  OF  TESTING 
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Vg  -  96.9  fps 

and  Che  on-board  inscrumencacion  must  be  shock  mounted  for  impact  velocities 
between  57  fps  and  96.9  fps  with  consideration  for  the  type  of  terrain  to  be 
encountered. 

CARGO  APPLICATION 

Example  2  illustrated  a  method  for  performing  high  altitude  parachute  tests 
at  lower  altitudes.  The  alternate  altitude  testing  proposition  may  also  be 
extended  to  heavy  cargo  weights  tested  at  low  altitude  and  velocity.  The 
possible  advantage  is  that  cargo  weights,  which  require  C-130  or  similar 
aircraft  as  launch  platforms,  can  be  reduced  to  sizes  which  can  be  flown  on 
aircraft  that  have  multiple  parent  station  wing  racks  utilizing  existing  test 
vehicles.  The  approach  to  testing  parachute  systems  at  higher  altitude  and 
velocity  with  lower  system  mass  is  similar  to  example  2,  but  in  reverse.  The 
method  of  analysis  is  similar  to  example  2  where  for  each  altitude  a  particular 
system  mass  and  test  velocity  are  determined  which  maintain  a  constant  BHR  and 
hence  maintain  a  constant  opening  shock  force  and  canopy  stress  distributions. 
Some  of  the  questions  to  be  addressed  are: 

1.  Does  Che  increase  in  test  altitude  introduce  unacceptable  complexity  in 
Che  daca-acquisicion  process? 

2.  How  is  the  required  test  altitude  affected  by  system  mass? 

3.  Is  the  required  test  dynamic  pressure  compatible  with  Che  available 
test  aircraft? 

Example  3  illustrates  a  typical  analysis. 

Example  3  -  Various  cargo  weights  of  4,000,  6,000,  8,000,  and  10,000  pounds  are 
Co  be  delivered  at  100  knots  near  sea  level  density  altitude.  The  retardation 
system  contains  a  solid  cloth  type  parachute  of  Do  135  feet  diameter. 

Determine  the  alternate  altitudes  for  aircraft  with  parent  rack  load  capability 
of  2,400  and  3,000  pounds.  The  parachute  geometry  is  taken  as: 


2a  b  b'  N 

—  “  0.668;  3  =  0.6214;  3  =  0.7806;  3  =  0.827 

D/,  a  a  a 


Therefore : 


a  *  45.090  ft;  b  =  28.019  ft;  b'  =  35.197  ft;  N  =  37.289  ft 

S  =  A  =  ^  =  14,314  ft^ 

o  so  4  o 
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CjjS^  «  0.75  X  »  10,735 


(12) 

(A14-31) 


k  -  1.46 
Cp^v  “ 

a.  Use  Equation  (A7-14)  to  determine  the  sea  level  inflation  distance. 


V  t 


8  O 


14W 


b.  Calculate  BHR  from  Equations  (A4-6) 

M  -  2W 

c.  Calculate  the  required  weight  and  inflation  distance  at  other 
density  altitudes. 


d.  Calculate  the  required  test  velocity  at  other  density  altitudes 
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e.  Calculate  the  refetence  inflation  time  at  the  particular  density 

altitude. 


'^ST 


£.  Calculate  the  time  of  occurrence  o f  the  maximum  shock  force  during 
the  inflation  process 


■(¥)’ 


g.  Calculate  the  maximum  shock  factor  during  Che  inflation  process 


X  . 

Lmax 


/21h\  ^ 

49  \  4  7 


h.  Calculate  the  maximum  shock  force 


V  ^ 

F  -  X.  r  -  X.  ^  SL  Cr^S 
max  imax  s  xroax  2  Do 

The  sea  level  characteristics  of  the  various  weights  are  tabulated  in  Table  5 
and  the  effects  of  density  altitudes  up  to  60,000  feet  are  plotted  in  Figures 
12,  13,  and  14. 

TABLE  5.  SEA  LEVEL  PARACHUTE  SYSTEM  CHARACTERISTICS  FOR  THE  VARIOUS  SYSTEM  WEIGHTS  OF 
EXAMPLE  3 


WEIGHT 

lb. 


INFLATION 

DISTANCE 

ft. 


16,962 

3,953 

2,018 


BMR 

t/to©x 

max 

5.748  X  lO-'* 

.43647 

3.698  X  1  r3 

.56943 

9.665  X  10-3 

.65311 

1.764  X  10-2 

.71187 

2.258  X  10-3 
1.113  X  10-3 
2.534  X  10-2 
4.249  X  10-2 


MAXIMUM 

FORCE 

lb. 


833 

4106 

9348 

15670 
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FIGURE  14.  EFFECT  OF  DENSTl/  ALTITUDE  ON  THE  REQUIRED  SYSTEM  TEST  VELOCITY  AND 

INFLATION  REFERENCE  TIME  FOR  A  CONSTANT  BALLISTIC  MASS  RATIO  FOR  EXAMPLE  3 


FIGURE  14.  EFFECT  OF  OENSTIY  ALTITUDE  OM  THE  REQUIRED  SYSTEM  TEST  VELOCITY  AND 

INFLATION  REFERENCE  TIME  FOR  A  CONSTANT  BALLISTIC  MASS  RATIO  FOR  EXAMPLE  3 
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The  weights  required  to  maintain  a  constant  BMR  for  altitudes  from  sea 
level  to  sixty  thousand  feet  are  presented  in  Figure  12.  Aircraft  with  parent 
station  racks  of  2,400  pounds  can  test  equivalent  sea  level  weights  of  4,000  and 
6,000  pounds  at  approximately  15,500  feet  and  26,500  feet,  respectively.  The 
sea  level  8,000-  and  10,000-pound  equivalents  require  altitudes  which  are  most 
likely  unacceptable.  An  increase  in  parent  rack  capability  to  3,000  pounds 
reduces  the  test  altitudes  for  the  sea  level  4,000-  and  6,000-pound  equivalent 
weights  to  approximately  9,000  and  21,000  feet,  respectively,  and  includes  the 
8,000-pound  weight  at  approximately  28,500  feet  test  altitude.  For  these 
altitudes,  long  focal  length  ground-to-air  photographic  coverage  and  one  or  more 
chase  aircraft  for  close  up  photographic  coverage  are  required.  The  question  of 
additional  cost  hinges  on  whether  the  photographic  equipment  and  chase  aircraft 
would  or  would  not  have  been  used  in  the  C-130  low  altitude  tests. 

The  effects  of  density  altitude  on  the  inflation  distance  are  presented  in 
Figure  13,  while  the  inflation  reference  time  and  test  velocity  are  presented  in 
Figure  14. 

The  problem  in  this  example  is  the  availability  of  an  aircraft  that  can  fly 
at  a  sea  level  dynamic  pressure  for  100  knots. 


IMPULSE  AND  MOMENTUM  DURING  PARACHUTE  INFLATION 

Another  approach  to  parachute  opening  dynamics  is  the  use  of  the  impulse 
imparted  to  the  system  by  the  varying  forces  during  the  time  of  canopy 
inflation.  In  this  section  a  method  has  been  developed  for  determining  the 
average  shock  factor  during  the  unfolding  and  elastic  stages  of  inflation.  The 
average  shock  factors  are  constant  with  altitude  for  constant  BMR  and  n.  The 
author  is  not  fully  satisfied  with  this  analysis  but  presents  the  method  as  a 
different  approach  to  the  use  of  canopy  inflation  impulse. 

For  deployments  where  the  BMR  is  lees  than  the  inflation  time,  tf,  is 
equal  to  the  reference  time,  t^.  When  the  BHR  exceeds  the  M^,  an  additional 
impulse  is  imparted  during  the  time  from  t^  to  t£  as  shown  in  Figure  15.  In 
the  general  case  where  n  is  not  zero,  the  impulse  in  the  interval  O^t^t^ 
is: 
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l-L  +  l.. 


TIME  t  (SEC) 


FIGURE  16.  IMPULSE  OF  THE  INFLATING  CANOPY  OF  EXAMPLE  2 


The  drag-area  and  velocity  ratios  may  be  expressed  as  functions  of  tin>e  ratio 
and  initial  drag-area  ratio.  Therefore  the  impulse  at  any  time  t^to  is  as 
follows: 

Using  for  the  general  case  in  the  unfolding  phase  of  inflation. 


-6  3 

+  2n(l-n)  +  n^l  dt 


+  2n(l-n) 


AC  any  time  ratio  the  average,  x^^y,  may  be  considered  as  the  area  under  the 
Xi  -  t  curve  divided  by  t„ 


NSWC  TR  85-24 


where 


(15) 


For  a  constant  BMR  and  n  with  altitude,  the  is  constant  for  all 

altitudes t  Once  the  x^^^  is  evaluated  at  the  initial  altitude  it  may  be  used 
at  any  other  altitude.  The  impulse  of  the  deployment  will  vary  since  t^ 
varies  with  altitude. 
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Whea  the  BMR  is  larger  than  an  additional  system  impulse  is  imparted 
between  Co^t<tf. 


f  (1)  “ 
y  lY.  j_  M'  jj  ‘  K' 

to  \^o  7M  L\to/ 


where 


and  Vq/V,  is  defined  by  equation  (18) 


At  t-t< 


1  /  (s)  Ao  vA 


■o  Vo  7M  IUo 


where 


‘  .  I  [(liY 

Vo  7M  L\to/  J 


The  time  ratio  tf/t^  is  also  constant  with  altitude 


^iav  Total  “  *iav 


^iav  Total 


<t/ld 
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^iav  Total 


(25) 


The  total  system  impulse,  shown  in  Figure  16,  may  be  considered  to  be  a 
rectangular  wave  form  which  is  more  convenient  to  use  as  a  basis  for 
calculations  than  the  usual  force-time  representation  of  the  inflating  parachute. 


IMPULSE  -  F,Xjg^ 


INFLATION  REFERENCE  TIME 
to  (SEC) 


FIGURE  16.  EFFECT  OF  ALTITUDE  AND  INFLATION  REFERENCE  TIME  ON  THE 
IMPULSE  OF  THE  INFLATING  CANOPY  OF  EXAMPLE  2 
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CONCLUSIONS 


The  theoretical  development  of  the  alternate  altitude  testing  technique  has 
demonstrated  the  feaseability  of  the  idea.  The  Ballistic  Mass  Ratio  meets  the 
criteria  for  a  genuine  performance  scale  factor  because  it  directly  or 
indirectly  affects  all  the  elements  which  define  the  parachute  inflation 
process.  Although  the  theory  was  directed  to  solid  cloth  parachutes,  it  also 
applies  to  other  decelerator  types.  For  other  decelerators  the  important  point 
is  to  express  the  parachute  inflation  distance  (V^Cq)  as  a  function  of 
density  altitude  and  adjust  the  required  weight  to  maintain  a  constant  BHR. 
Testing  at  a  constant  dynamic  pressure  and  a  constant  BMR  at  any  density 
altitude  produces  the  same  maximum  opening  shock  force  and  stress  distribution. 
The  time  ratio  (t/tg)  at  the  time  of  the  maximum  shock  force  occurence  during 
canopy  inflation  remains  constant,  although  the  parachute  inflates  more  rapidly 
as  altitude  increases.  When  viewed  from  the  time  ratio  aspect,  the  system 
geometry,  instantaneous  shock  factor,  velocity  decay,  dynamic  pressure 
variation,  gore  pressure  distribution,  average  pressure  coefficients  and  stress 
distribution  are  the  same  for  all  altitudes. 

Limitations  to  alternate  altitude  testing  are  evident  in  the  extreme  weight 
and  low  velocity  required  in  example  2  and  the  low  dynamic  pressure  required  in 
example  3.  However,  there  is  a  wide  range  of  applications  between  these 
examples  where  alternate  altitude  testing  should  be  a  viable  approach. 

A  method  of  addressing  the  impulse  and  momentum  of  an  inflating  parachute 
canopy  was  also  derived.  The  average  shock  factor  during  canopy  inflation  was 
shown  to  be  a  function  of  the  ballistic  mass  ratio. 
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A  TCCHKIQUE  FOR  THE  CALCULATION  C?  THE  CPENII.'C.-SHOCK  FORCES  FOR 

SHR.^''R;-L  types  of  solid  cloth  p.-eachutes 

W.  P.  L'jd*ke 

I.'aval  Ordr.'ir.ce  Laboratory 
Sliver  Sprlr.g,  Karylar.d 


Abstract 

An  analytical  method  of  calculating 
parachute  opening-shock  forces  based  upon 
vjind-tunnel  derived  drag  area  time  signa¬ 
tures  of  several  solid  cloth  parachute 
tj'pes  in  conjunction  with  a  scale  factor 
and  retardation  system  steady-state  param¬ 
eters  has  been  developed.  Methods  of 
analyzing  the  inflation  tlm.e,  geometry, 
cloth  airflow  properties  and  materials 
elasticity  are  Included.  The  effects  of 
mass  ratio  and  altitude  on  the  magnitude 
and  time  of  occurrence  of  the  maximum 
opening  shock  are  consistent  with  observed 
field  test  phenomena, 

I.  Introduction 


altitude  increases,  effects  of  altitude  or, 
opening-shock  force,  finite  a-d  infinilc 
mass  operation,  and  inflation  distance. 

II.  Development  of  Velocity  Ratio  and 
Force  Ratio  Equations  During 
the  Unfolding  Phase  of 
_ Parachute  Deployment  ■ 

The  parachute  deployment  would  take 
place  Ir.  a  horizontal  attitude  in  accord¬ 
ance  with.  Newton '  8  '  second  lav;  of  motion. 

DP  =  ma 


dt 


In  1965 »  the  Naval  Ordnance  Laboratory 
(NOL)  v.'as  engaged  in  a  pro.iect  which 
utilized  a  35-foot-diameter,  10-percent 
extended-skirt  parachute  (type  T-10)  as 
the  second  stage  of  a  retaraation  system 
for  a  250-pound  payload.  Deployment  of 
the  T-10  parachute  was  to  be  accomplished 
at  an  altitude  of  100,000  feet.  In  this 
rarefied  atmosphere,  the  proble.m  v;as  to 
determine  the  second  stags  deployment 
conditions  for  successful  operation.  A 
search  of  available  field  test  information 
indicated  a  lack  of  data  on  the  use  of 
solid  cloth  parachutes  at  altitudes  above 
30,000  feet. 

The  approach  to  this  problem  was  as 
follows:  Utilizing  existing  v;lr.d-tunnel 
data,  lov;-altltude  field  test  data,  and 
reasonable  assumptions,  a  unique  engi¬ 
neering  approach  to  the  inflation  time 
and  opening-shock  problem  was  evolved  that 
provided  satisfactory  results.  Basically, 
the  method  combines  a  wind-tunnel  derived 
drag  area  ratio  signature  as  a  function 
of  deployment  time  with  a  scale  factor  and 
Newton's  second  law  of  motion  to  analyze 
the  velocity  and  force  profiles  during 
deployment.  The  parachute  deployment 
sequence  is  divided  into  tvro  phases.  The 
first  phase,  called  "unfolding  phase," 
vrtiere  the  canopy  is  undergoing  changes  in 
shape,  is  considered  to  be  Inelastic  as 
the  parachute  Inflates  initially  to  its 
steady-state  aerodynamic  size  for  the 
first  time.  At  this  point,  the  "elastic 
phase"  is  entered  v;here  it  is  co:-.sldered 
that  the  elasticity  of  the  parachute 
materials  enters  the  problem  and  resists 
the  applied  forces  until  the  canopy  has 
reachied  full  Inflation. 

The  developed  equations  are  in  agree¬ 
ment  with  the  observed  performance  of 
solid  cloth  parachutes  In  the  field,  such 
as  the  decrease  of  inflation  time  as 


It  was  recognized  that  other  factors,  such 
as  included  air,  mass,  apparent  mass,  and 
their  derivatives,  also  contribute  forces 
actl'  1  on  the  system.  Since  definition 
of  ti.ese  parameters  was  difficult,  the 
analysis  was  conducted  in  the  simplified 
form  showTi  above.  Comparison  of  calculated 
results  and  test  results  indicated  that 
the  omitted  terms  have  a  small,  effect. 


t  V 


.Multiplying  the  right-hand  side  of  equa¬ 
tion  (1)  by 


1  = 


s  o  D  o 

V„t,^CnS„ 
S  0  D  0 


and  rearrangihg 


"o  Cd^o 


dt 


-2v; 


Pgv  t  C-S 
so  Do 


(2) 


In  order  to  integrate  the  left-hand  term 
of  equation  (2),  the  drag  area  ratio  must 
be  defined  for  the  type  of' parachute  under 


A-2 


i;'.E.lysia  aa  a  i'unetlon  of  deployment 
time,  tjj. 


FIG.  1  TYPICAL  INFINITE  MASS  FORCE-TIME  HISTORY  OF  A 
SOLID  CLOTH  PARACHUTE  IN  A  WIND  TUNNEL 
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REPRODUCED  FROM  REFERENCE  (1) 


FIG.  3  TYPICAL  FORCE-TIME  CURVE  FOR  A  10%  EXTENDED 

SKIRT  PARACHUTE  UNDER  INFINITE  MASS  CONDITIONS. 


FIG.  4  TYPICAL  FORCE-TIME  CURVE  FOR  A  PERSONNEL 
GUIDE  SURFACE  PARACHUTE  UNDER  INFINITE 
MASS  CONDITIONS 


Figure  1  illustrates  a  typical  solid 
cloth  parachute  v;ind-tunnel  infinite  mass 
-force-time  history  after  snatch.  In 
infinite  mass  deployment,  the  maximum  size 
and  maximum  shock  force  occur  at  the  time 
of  full,  inflation,  .  tf.  ,  Hov(eyer,  .  tf  is 
inepproprlate  -for;  analyslPsincei  it  :1s'  ’ ' 
dependent  .  upon-,  the  applied  load; ' ;  structural 
strength,  and  materials  elasticity. 4.  The  ? 
reference-  time,  tQ,  ;|v;here  -the  -parachute ;?has 
attained  its  steady-state,  aerodynamic  size ; 
for  the  first  time,  is .used  as  the  basis  for 
performance  calculations. 


At  any  instant  during  the  unfolding 
phase,  the  force  ratio  F/Pg  can  be  deter¬ 
mined  as  a  function  of  the  time  ratio, 
t/ f  c  . 


S  .,  3  Do 


Since  the  wind-tunnel  velocity  and  density 
are  constant  during  infinite  mass  deployment 

F  _ 


0  ’  a04  0.08  0.12  0.18 

TIME  (SEC) 

FIG.  5  TYPICAL  FORCE-TIME  SIGNATURE  FOR  THE  ELLIPTICAL 
PARACHUTE  UNDER  INFINITE  MASS  CONDITIONS 


TIME  (SEC) 


FIG,  2  TYPICAL  FORCE-TIME  CURVE  FOR  A  SOLID  FLAT  PARA¬ 
CHUTE  UNDER  INFINITE  MASS  CONDITIONS. 


FIG.  6  TYPICAL  FORCE-TIME  SIGNATURE  FOR  THE  RING 

SLOT  PARACHUTE  20%  GEOMETRIC  POROSITY  UNDER 
INFINITE  MASS  CONDITIONS 


In-'ir.tte  mass  opening-shock  signatures 
cf  several  types  of  parachutes  are  pre¬ 
sented  in  Figures  2  through  6.  Analysi.s 
of  these  signatures  using  the  force  ratio, 
r/Fg,  -  ti.r.e  ratio,  t/t^,  technique 
indicated  a  similarity  In  the  performance 
of  the  various  solid  cloth  types  of 
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parachutes  which  v;ere  examlr.ed .  Hie  geo¬ 
metrically  porous  ring  slot  parachute 
displayed  a  completely  different  signature, 
as  v;a8  expected.  These  data  are  illus¬ 
trated  in  Figure  ?•  If  an  initial  boundary/ 
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FIG.  7  DRAG  AREA  RATIO  VS.  TIME  RATIO 


condition  of  Cj^S/C^S^  =  0  at  time  t/t^  =  0 

Is  assu.T.od,  then,  the  data  can  be  approxi¬ 
mated  by  fitting  a  curve  of  the  form 


c  c 

c.^s 

Do 


(3) 


At  the  time-  that  equation  (5)  i--'as  ascer¬ 
tained,  It  suggested  that  the  geometry  of 
the  deploying  parachute  •.•;as  independent 
of  density  and  velocity.  It  v.-as  also 
postulated  that  although  this  expression 
had  been  determined  .'ci'  the  infinite  mass 
condition,  it  would  also  be  true  for 
the  finite  mass  case.  'Ihis  phenomenon 
has  since  been  independently  observed  and 
confirmed  by  Berndt  and  Be  ’..'eese  in 
reference  {it). 

Since  the  drag  area  ratio  was  determined 
from  actual  parachute  deployments,  it  was 
assu-med  that  the  effects  of  apparent  mass 
and  included  mass  on  the  deploi-ment  force 
history  were  accommodated. 

The  right-hand  term  of  equation  (2) 
contains  the  expression 


2W 


(6) 


This  term  can  be  visualised  as  shoivn  in 
Figure  8  to  be  a  ratio  of  the  retarded 
mass  (including  the  parachute)  to  an 
associated  mass  of  atmosphere  contained 
in  a  right  circular  cylinder  which  is 
generated  by  m’sving  an  inflated  parachute 
of  area  CjSq  ^  distance  equal  to  the 
product  of  Vgto  through  an  atmosphere  of 
density,  o. 


A  more  realistic  drag  area  ratio  expres¬ 
sion  was  determined  which  Includes  the 
effect  of  initial  area  at  line  stretch. 


FIG.  8  VISUALIZATION  OF  THE  MASS  RATIO  CONCEPT 


The  mass  ratio,  M,  Is  the  scale  factor 
which  controls  the  velocity  and  force 
profiles  during  parachute  deploiTnent. 
Substituting  M  and  Cj^S/CqS^  into  equa¬ 
tion  (2),  integrating,  and  solving  for 


vfhero  n  is  the  ratio  of  the  projected 
mouth  area  at  line  stretch  to  the  steady- 
state  projected  frontal  area.  Expanding 
equ.'ition  (4) 


V. 


1 


(1 


^  \% 


nd  -T?)(  t 


(7) 
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The  Instantaneous  shock  factor  Is  defined 
as 


F_  1 


«  l^^s  Vo 


Vo 


(r) 


from  equations  (5)  and  (7) 


r,  <8) 


III.  Maximum  Shock  Force  and  Time  of 
Occurrence  During  the  Unfolding  Phase 


The  time  of  occurrence  of  the  maximum 
Instantaneous  shock  factor,  x^.  Is  diffi¬ 
cult  to  determine  for  the  general  case. 
However,  for  n  =  0,  the  maximum  shock 
factor  and  time  of  occurrence  are  readily 
calculated.  For  n  =■  0 


, .  (i)^  ^ 


Setting  the  derivative  of  x^^  with  respect 

to  time  equal  to  zero  and  solving  for 

t/t-  at  X. 
o  1  max 


f-s-'  .  (mV 

\to/  @  Xj_  \  ^  / 


and  the  maximum  shock  factor  Is 


i  max 


U 

16  /21m\ 

^9  U  / 


Equations  (9)  (10)  are  valid  for  value., 

of  H  s  ~  (0.19),  since  for  larger  values 

of  M,  the  maximum  shock  force  occurs  In 
the  elastic  phase  of  inflation. 


Figures  9  and  10  illustrate  the  velocity 
and  force  profiles  .generated  from  equat.lon.s 
(7)  and  (8)  for  initial  projected  area 
ratios  of  n  =  0,  and  0.2  with  various  m.ass 
ratios. 


Methods  for  Calculation  of 
the  Reference  Time,  tp 


If  the  altitude  variation  during  deployment 
is  small,  then,  the  density  may  be  consid¬ 
ered  as  constant 


The  ratio  concept  is  an  ideal  method  to 
analyze  the  effects  of  the  various  param¬ 
eters  on  the  velocity  and  force  profiler, 
of  the  opening  parachutes;  hovievcr,  a 
means  of  calculating  t^is  required  before 
specific  values  can  be  co.mputed.  Methods 
for  computing  the  varying  mass  flov,-  into 
the  inflating  canopy  mouth,  the  varying 
mass  flow  out  through  the  varying  inflated 
canopy  surface  area,  and  the  volume  of 
air,  which  must  be  collected  during 
the  inflation  process  are  required. 


Figure  11  represents  a  solid  cloth- 
type  parachute  canopy  at  some  instant 
during  inflation.  At  any  given  instant, 
the  parachute  drag  area  is  proportional 
to  the  maximum  Inflated  diameter.  A].so, 
the  maximum  diameter  in  conjunction  with 
the  suspension  lines  determines  the  inflow 
mouth  area  (A-A)  and  the  pressurized 
canopy  area  (B-B-B).  This  observation 
provided  the  basis  for  the  following 
assximptions.  The  actual  canopy  shape  is 
cf  minor  importance. 


a.  The  ratio  of  the  instantaneous 
mouth  inlet  area  to  the  steady-state 
mouth  area  is  In  the  same  ratio  as  the 
instantaneous  drag  area. 


b.  The  ratio  of  the  instantaneous 
pressurized  cloth  surface  area  to  the 
canopy  surface  area  is  in  the  sa.Tie  ratio 
as  the  instantaneous  drag  area. 


c.  Since  the  suspension  lines  in  the 
ur.pressurized  area  of  the  canopy  are 
straight,  a  pressure  differential  has  not 
developed,  and,  therefore,  the  net  air¬ 
flow  in  this  zone  is  zero. 


Based  on  the  foregoing  assumptions, 
the  mass  flov.’  equation  can  be  v/ritten 


dm  =  m  Inflovf  -  m  out  flov,' 


p  -^  =  pVA,.  -  pA  P 
dt  ‘‘  ® 


V/V  VELOCITY  RATIO  V/V.  VELOaTY  RATIO 
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dV 


CpS 


CrS 

D 


p  —  =  PVA^  — -  pf,  — - 
dt  °  “O  z’  e 


Cd^o 


(11) 


FIO.  11  PARTIALLY  INFLATED  PARACHUTE  CANOFV 


From  equation  (3) 


C„S 
C 


D^o  Vo/ 


i  for  n  “  0 


From  equation  (7) 


1  + 


-f-T  ’ 


;  Ti  =  0 


From  equation  (26) 


P  k 


„?n 


-Aso*'- 


Pn 

dt 

(1^) 


Integrating: 


^o 


-Aso'' 


+ 


2n 

dt 

(13) 


Measured  values  of  n  Indicate  a  data 
range  from  0.57^  through  £>.771.  A 
convenient  solution  to  the  reference  time 
equation  evolves  when  n  is  assigned  a 
value  of  1/2,  Integrating  equation  (13) 
and  using 


Vat„M 

s  o 


2W 


gpCpS^ 


LET  =  - 


6PZo 


2W 


Cr«S 

Do 


^«Io  ■ 


iW- 


^o  " 


14W 


epVgC 


D''o 


Kl 


-1 


{lU) 


Equal lor.  (l4)  expresses  the  unfolding 
reference  time,  to.  In  terms  of  mass, 
altitude,  snatch  velocity,  airflow  char¬ 
acteristics  of  the  cloth,  and  the  steady- 
state  parachute  geometry.  Mote  that  the 
term  g^V^j/W  Is  the  ratio  of  the  Included 
air  mass  to  the  mass  of  the  retarded 
hardware.  Multiplying  both  sides  of 
equation  (14)  by  Vg  demonstrates  that 


V  t  •»  a  constant  v/hlch  Is 

a  function  of  altitude 


Figures  12  and  13  Indicate  the  para¬ 
chute  unfolding  time  and  unfolding 
distance  for  values  of  n  *•  1/2  arid  n  = 
0.63246,  Mote  th(.  .  'Elation  and  conver¬ 
gence  v/lth  rising  altitude.  The  opening- 
shock  force  is  strongly  influnncnd  by 
the  Inflation  time.  Because  of  this,  the 
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value  of  Iq  calculated  by  using  a  real¬ 
istic  value  of  n  should  be  used  in  the 
lower  atmosphere. 

As  or.  example  of  this  method  of  opening- 
shock  analysis,  let  us  examine  the  effect 
of  altitude  on  the  opening-chock  force  of 
a  T-lO-type  parachute  retarding  a  POO- 
pourid  weight  from  a  snatch  velocity  cf 
Vg  -  IJOO  feet  per  second  at  sea  level. 
Conditions  of  constant  velocity  and  con¬ 
stant  dynamic  pressure  are  investigated. 

The  results  are  presented  in  Figure  l4. 

At  low  altitudes,  the  opening-shock  force 
is  less  than  the  steady-state  drag  force; 
however,  as  altitude  rises,  the  opening 
shock  eventually  exceeds  the  steady-state 
drag  force  at  some  altitude.  This  trend 
is  in  agreement  with  field  test  observatlors. 

V.  Correction  of  t^,  for 
Initial  Area  Effects 

Ihe  unfolding  reference  time,  t,,, 
calculated  by  the  previous  methods  assumes 
that  the  parachute  inflates  from  sero  drag 
area.  In  reality,  a  parachute  has  a  drag 


CLOTH  MIL  C  7020,  TYPE  III 

3S  FOOT;  10*/.  EXTENDED  SKIRT  CANOPY 


UNFOLOINO  DISTANCE  IPEETI 


no.  1)  EFFECT  OF  ALTITUDE  ON  THE  UNFOLDING  DIS¬ 
TANCE  AT  CONSTANT  VELOCITY  AND  CONSTANI 
DYNAMIC  PREnURf  FOR  n  •  V2  AND  n  •  O.S124a 


100  200 


FORCE  (THOUSANOSOF  FOUNDS) 

FlO.  14  VARIATION  OF  STEADY-STATE  DRAG.  F^  AND  MAXI¬ 
MUM  OPENINO  SHOCK  WITH  ALTITUDE  FOR  CONSTANT 
velocity  and  constant  dynamic  PRItSURS 


UNFOLOINO  TIME  <SEC.I 


FIG.  12  EFFECT  OF  ALTITUDE  ON  THE  UNFOLOINO  TIME 
AT  CONSTANT  VELOCITY  ANO  CONSTANT 
DYNAMIC  PRESSURE  FOR  n  •  1/2  ANO  n  •  0  62206 


aiea  at  the  beginning  of  inflation.  Once 
to  has  been  calculated,  a  correction  can 
be  applied,  based  upon  what  is  known  about 
the  initial  conditions. 
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C*ico  A  -  the  initial  projected 

area  Is  kno’-rn 


than  the  maximum  extensibility,  c  ,  of 
the  materials. 


t 


1 


^calculated 


Development  of  the  analysis  In  the 
elastic  phase  of  inflation  is  similar  to 
the  technique  used  in  the  unfolding  phase, 
ilewtor. *8  second  law  of  motion  is  used, 
together  with  the  drag  area  ratio  signa¬ 
ture  and  inass  ratio 


corrected 


^calculated 


(15) 


Case  B  -  V/her,  the  Initial  drag  area 
Is  known 

Vo  Vo/ 


(16) 

The  mass  ratio  should  row  be  adjusted  for 
the  corrected  before  velocity  and  force 
profiles  are  de.erp.lr.ed. 

VI.  Opening-Shock  Force,  Velocity  Ratio, 
and  Inflation  Time  During  the 
Elastic  Phase  of  Parachute  Inflation  . 


The  mass  ratio,  M,  Is  an  Important 
parameter  In  parachute  analysis.  For 
values  of  M  <<4/21,  the  maximum  opening- 
shock  force  occurs  early  In  the  inflation 
process,  and  the  elastic  properties  of 
the  canopy  are  rot  significant.  As  the 
mass  ratio  approaches  M  -  4/21,  the 
magnitude  of  the  opening-shock  force 
Increases,  a.-.d  the  time  of  occurrence 
happens  later  in  the  deployment  sequence. 
For  mass  ratios  M  >  4/21,  the  maximum 
shock  forC'C  v;lll  occur  after  the  reference 
time,  to.  Parachutes  designed  for  high 
mass  ratio  operation  must  provide  a 
structure  of  sufficient  constructed 
strength,  Fp,  so  that  the  actual  elon¬ 
gation  of  trie  canopy  under  load  Is  less 


A-hlch  Is  still  valid,  as  shown  In  Figure  7 


V2 


Integrating  and  solving  for  — 

Vs 


vfhere  —  is  the  velocity  ratio  of  the 

s 

u-nfoldlng  process  at  time  t  •  t^. 


1  + 


1 

M 


2 


The  Instantaneous  shock  factor  In  the 
elastic  phase  becomes 


The  end  point  of  the  Inflation  process 
depends  upon  the  applied  loads,  elasticity 
of  the  canopy,  and  the  constructed  strength 
of  the  parachute.  A  linear  load  elongation 
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relationship  Is  utilized  to  determine  the 
maximum  drag  area. 


e  *roax 


F 


F 


c 


Tne  r.ext  force  In  the  scries  at  conctant  q 


F 


1  = 


5^„F.  -t 


c 


F  c 


max 


where 


(20) 


(1  + 


Ihc  force,  F,  is  initially  the  instanta¬ 
neous  force  at  the  end  of  the  unfolding 
process 


O  8 


(21) 


Subsequent  elongations  ir.  the  system  can 
be  shov.T.  to  be 

*1  “  'o  +  'o)^ 


2  2 

where  X©  is  the  shock  factor  of  the  *2  “  *o  *o  ®o^  ^ 

unfolding  phase  at  t  »  t^ 


1^1  + 1  [g.yal!  +  na.r  p)  +  ^2]J 

(22) 

Since  the  inflated  shape  Is  defined,  the 
drag  coefficient  is  considered  to  be 
constant,  and  the  Instantaneous  force  is 
proportional  to  the  dynamic  pressure  and 
projected  area.  The  maximum  projected 
area  would  be  developed  if  the  dynamic 
pressure  remained  constant  during  the 
elastic  phase.  Under  very  high  mass 
ratios,  this  is  nearly  the  case  over  this 
very  brief  time  psriod;  but  as  the  mass 
ratio  decreases,  the  velocity  decay  has 
a  more  significant  effect.  The  simplest 
approach  for  all  mass  ratios  is  to 
determine  the  maximum  drag  area  of  the 
canopy  as  if  elastic  inflation  had 
occurred  at  constant  dynamic  pressure. 

Then  utilizing  the  time  ratio  determined 
as  an  end  point,  intermediate  shock 
factors  can  be  calculated  from  equation 
(19)  and  maximum  force  assessed. 

The  initial  force,  X^Fg,  causes  the 
canopy  to  increase  in  projected  area.  The 
new  projected  area  in  turn  increases  the 
total  force  on  the  canopy  which  produces 
a  secondary  projected  area  Increase.  The 
resulting  series  of  events  are  resisted 
by  the  parachute  materials.  The  parachute 
must,  therefore,  be  constructed  of  suffl- 
clci.t  strength  to  prevent  the  elo.'igatlon 
of  the  materials  from  exceeding  the 
maximum  elongation. 


*0  “  p  "  *max  (23) 


The  required  canopy  constructed  strength 
can  be  determined  for  a  given  set  of 
deployment  co.nditlons.  Tne  limiting  value 
of  the  series  (c^)  determines  the  end 
point  time  ratio. 


^U^max 


(1 


‘t) 


2 


(1  +  «/.) 


1/3 


(24) 


Figure  13  illustrates  the  maximum  drag 
area  ratio  as  a  function  of  e^. 


ftc.  IS  MAXIMUM  DRAG  AREA  RATIO  VS  INITIAL 
ELONGATION 


A-10 


NSWC  TR  85-24 


VII,  Applicatio:'.  of  Cloth  Permeability 

to  the  Calculatio".  of  the  Inflation 
Tlr^e  of  Solid  Cloth  Parachutes _ 

The  mass  outflov;  through  the  pressur¬ 
ized  reglOT;  of  an  inflating  solid  clot-i 
parachute  at  any  instant  is  dependent  upon 
the  canopy  area  which  is  subjected  to  air- 
flow  and  the  rate  of  airflow  through  that 
area.  The  variation  of  pressurized  area  os 
a  function  of  reference  time,  t^,  t;ar. 
earlier  asstuned  to  be  proportional  to  the 
Instantaneous  drag  area  ratio,  leaving  the 
rate-of-airflo'.-f  problem  to  solve.  Tnc 
permeability  parameter  of  cloth  was  a 
natural  choice  for  dotenrlning  the  rate  of 
airflow  through  the  cloth  as  a  function  of 
pressure  differential  across  the  cloth. 
Heretofore,  these  data  have  been  more  of 
a  qualitative,  rather  than  quantitative, 
value.  A  new  method  of  analysis  was 
developed  './herein  a  generalized  curve  of 
the  fonr,  p  =  k(dP)'^  v/as  fitted  to  cloth 
permeability  data  for  a  number  of  different 
cloths  and  gives  surprisingly  good  agree¬ 
ment  over  the  pressure  differential  range 
of  available  data.  The  pressure  differ¬ 
ential  was  then  related  to  the  trajectory 
conditions  to  give  a  generalized  exprcsslo.-. 
v/hlch  can  be  used  in  the  finite  mass  ratio 
range,  as  v;ell  as  the  infinite  mass  case. 
The  permeability  properties  were  trans¬ 
formed  into  a  nass  flow  ratio,  M*,  which 
shows  agreement  v/lth  the  effective  porosity 
concept. 

Measured  and  calculated  permeability 
pressure  data  for  several  standard  cloths 
are  illustrated  in  Figure  l6.  This  method 
has  been  applied  to  various  types  of  cloth 
between  the  extremes  of  a  highly  permeable 
3-nonsne  silk  to  a  relatively  impervious 
parachute  pack  container  cloth  with 
reasonably  good  results,  see  Figure  17. 

The  ca-nopy  pressure  coefficient, 
is  defined  as  the  ratio  of  the  pressure 
differential  across  the  cloth  to  the 
dynamic  pressure  of  the  free  stream. 


C  =  ^  o  PC-rs-err-al )  -■  P(exter.nal) 

**  q  1/2  pV^ 

(23) 

v/here  V  is  based  on  equation  (7). 

The  pzrmeablllty  expression,  P  =  k(aP)'^ 
becomes 


(2G) 


Altho'jgh  some  progress  has  been  made 
by  .Melzig  and  others  or  the  measurement  of 
the  variation  of  the  pressure  coefficient 
on  an  actual  inflating  canopy,  this 
dimension  a’'d  its  variation  v/lth  time  are 
still  dark  areas  at  the  time  of  this 
writing.  At  the  present  time,  a  constant 
average  value  of  pressure  coefficient  is 


used  in  these  calculations.  Figure  ifl 
presents  the  effect  of  pressure  coefficient 
and  altitude  or,  the  ur:folding  tine  for 
constant  deployment  conditioriS. 

It  is  well  known  that  the  ii.fl.atlon 
time  of  solid  cloth  parachutes  decreases 
as  the  operatlo.nal  altitude  Increases. 

Tnis  effect  can  be  explained  by  consid¬ 
ering  the  ratio  of  the  mass  outflow 
through  a  unit  cloth  area  to  the  mass 
inflov/  through  a  unit  mouth  area. 

K*  =  mass  ilo./  ratio  - - 

macs  inflow 

where 

mass  outflow  P.  -5^°-(per  cloth  urea) 
“  ft'^-sec 


and 

ina.S8  Inflow  =  V  ^>^2  area) 

P  ft^-sec 


REPROOUCCO  FROM  REFERENCE  (4) 


Fic.  IS  nominal  porositv  of  rarachutf  material 

VS  DIFFERENTIAL  PRESSURE. 
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Therefore,  the  nase  flov/  ratio  becomes 


AVERAGE  CANOPY  RhESSURE  COEFPICIENT 
DURING  INFLATION  INCLUDING  THE  VENT 


VD  V 


"n".  1/2 
Ea  (14) 


REALISTIC  •  n' 
EO  <l» 


-  -  Kt) 


Effective  porosity,  C,  la  defined  as 
the  ratio  of  the  velocity  through  the 
cloth,  u,  to  a  fictitious  theoretical 
velocity,  V,  v/hlch  v/111  produce  the  par¬ 
ticular  IPs  l/Pov*-. 


o  ask 


o .  ft* 

*  -* 

5  ’ 

^  as 
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60  a  FT 


JOILII- 
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FIG.  It  EFFECT  OF  FREtSURi  COEFFICIENT  AND  ALTITUDE 
ON  THE  UNFOLDINO  TIME. 
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effective  porosity,  C  — 

V 


Comparison  of  the  mass  flow  ratio  and 
previously  published  effective  porosity 
data  Is  shown  in  Figure  19.  Tne  effecta 
of  altitude  and  velocity  on  the  mass  flow 
ratio  are  presented  In  Figures  20,  and  21 
for  constant  velocity  and  constant 
altitude.  Tne  decrease  of  cloth  perme¬ 
ability  with  altitude  Is  evident. 

TJie  permeability  constants  "k"  and  "n" 
cam  be  determined  from  the  permeability 
pressure  differential  data  as  obtained 
from  an  instrument  such  as  a  Frazier 
Permeameter.  Two  data  points,  "A"  and 
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"B, "  are  selected  ir.  such  a  maimer  that 
point  "a"  is  in  a  lov.-pressure  zone  belov; 
the  knee  of  the  curve,  and  point  "Ft"  is 
located  ir  the  upper  er.d  of  the  high- 
pressure  cone,  as  s'r.o..T.  in  Figure  f2. 

The  two  standard  '’.easurements  of  1/? 
inch  of  water  and  ?0  inches  of  water 
appear  to  be  good  data  points  if  both  are 


riG.  22  LOCATION  OF  DATA  POINTS  FOA  DETERMINATION 
OF  "k~  AND"n" 


available  on  the  same  sample.  Substi¬ 
tuting  the  data  from  points  "A"  and  B" 

into  P  =  k(AP)''  : 


FIG.  20  EFFECT  OF  ALTITUDE  ON  MASS  FLOW  RATIO  AT 
CONSTANT  VELOCITY 


FIG  21  EFFECT  OF  VELOCITY  ON  MASS  FLOW  RATIO  AT 
CONSTANT  DENSITY 


Vill.  Determination  of  the 
Parachute  Included  Volume 
and  Associated  Air  Mass 


Before  the  reference  time,  to,  and 
inflation  time,  tf,  can  be  calculated, 
the  volume  of  atmosphere,  Vo ,  which  Is 
to  be  collected  during  the  Inflation 
process  must  be  accurately  known.  This 
requirement  dictates  that  a  realistic 
Inflated  canopy  shape  and  associated 
volume  of  atmosphere  be  determined. 

Figure  ?■’  reproduced  from  reference 
(5)'  The  technique  of  using  lampblack 
coated  plates  to  determine  the  airflow 
patterns  around  metal  models  of  inflated 
canopy  shapes  was  used  by  the  investigator 
of  refere.nce  (5)  to  study  the  stability 
c'r.arac Toristlcs  of  contemporary  parachutes, 
l.e.,  19^3-  A  by-product  of  this  study 
is  that  it  is  clearly  shown  that  the 
volume  of  air  within  the  canopy  bulges 
out  of  the  ca.nopy  mouth  (indicated  by 
arrows)  and  extends  ahcaa  of  the  canopy 
hem.  This  volume  must  be  collected  durljig 
the  Inflation  process.  Another  neglected, 
but  significant,  source  of  canopy  volume 
exists  in  the  billowed  portion  of  the  gore 
panels , 
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VENT  PARACHUTE 


REPRODUCED  FROM  REFERENCE  (S) 

FIG.  23  AIRFLOW  PATTERNS  SHOWING  AIR  VOLUME 
AHEAD  OF  CANOPY  HEM 


The  steady-state  canopy  shape  has  been 
observed  In  wind-tunnel  and  field  tests  to 
be  elliptical  in  profile.  Studies  of  the 
Inflated  shape  and  included  volume  of 
several  parachute  types  (flat  circular, 

10  percent  extended  skirt,  elliptical, 
hemispherical,  ring  slot,  ribbon,  and 
cross)  are  documented  in  references  (6) 
and  (7).  These  studies  demonstrated  that 
the  steady- state  profile  shape  of  Inflated 
canopies  of  the  various  types  can  be 
approximated  to  be  two  ellipses  of  common 
major  diameter,  2a,  and  dissimilar  minor 
diameters,  b  and  b',  as  shown  in  Figure  24. 
It  was  also  shown  that  the  volume  of  the 
ellipsoid  of  revolution  formed  by  revolving 
the  profile  shape  about  the  canopy  axis 
VUG  a  good  approximation  of  the  volume  of 
atmosphere  to  be  collected  during  canopy 
inflat..on  and  Included  the  air  volume 
extend  cl  ahead  of  the  parachute  skirt  hem 
togethei  v;ith  the  billovfed  gore  volume. 


FIG.  24  PARACHUTE  CROSS  SECTION  NOMENCLATURE 


Tables  I  and  II  are  summaries  of  tost 

results  reproduced  from  references  (6)  a.nd 

(7),  respectively,  for  the  convenience  of 

the  reader. 
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TABLE  I  SUMMARY  OF  PARACHUTE  SHAPE  TEST  RESULTS 
FOR  12  GORE  AMD  16-60RE  CONFIGURATIONS 


Paraefiut* 

TVP* 

No.  o1 

Gortt 

Soiiwntion 
Lnm  LMiffh 
•nehtt 

ValociVY 

mp(t  fpt 

2i 

Oo 

Scil«  Pactof.  K 
2i  2a 

0. 

2a 

L 

N 

a 

Aatt  Rilto 

b  b’ 

a  a 

b  b’ 
i  *: 

VoHima  in 

Vm  Vj 

3 

''m 

■  “  ■ 

12 

34 

60 

73 

.645 

690 

656 

6115 

8617 

1  4932 

4476 

4461 

6960 

156 

16 

34 

60 

73 

663 

666 

620 

.5556 

.9039 

1.459/ 

445C 

4100 

7325 

1  55 

10%  Cattndtd 

12 

34 

too 

147 

663 

.652 

661 

6424 

8660 

19214 

3926 

4400 

6783 

173 

$h«i 

16 

34 

17 

25 

664 

.640 

785 

5560 

8502 

1.4062 

4051 

3920 

6197 

1  53 

12 

KB 

■a 

916 

612 

5626 

9667 

1.5263 

3322 

5606 

16 

■a 

la 

675 
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.6169 

6163 

1  4332 

2726 

4405 

HomitphtfiCf' 

12 

n 

M 

996 

1  254 

1  0005 
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1  9065 

6224 

8666 

16 

la 

Ba 

994 

1  165 

9129 

9360 

1.6509 

5521 

6370 

Rtn^oi 

12 

34 

26 

37 

607 

654 

653 

.6566 

1735 

1  530 
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16%  GoomaOic 

12 

34 

too 

147 

616 

663 

.922 

6566 
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16 

25 

37 
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62? 

6004 

6690 

1  4694 

3600 
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5565 

1  50 

16 
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664 

6004 

•690 

1  4694 

3600 

3965 

6030 
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16 

KB 
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644 

6004 

•690 

1  4694 

3600 

4430 

660^ 

1  62 

Ribboo 

12 

34 

25 

37 

566 

632 

.659 

6556 

6766 

1.9326 

3600 

3323 

SJ» 

1  40 

24%  Q«om«trM 

12 

34 

100 

147 

615 

663 

.637 

65S6 

8766 

1  5326 

3«00 

3714 

6163 

1  52 

Porowty 

12 

34 

200 

293 

632 

661 
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6556 

8766 

1  5326 

3600 

4260 

6683 

1.76 

16 

34 

25 

37 

603 

660 

79? 

5570 

8578 

1.4146 

MM 

3436 

5356 

1  41 

16 

34 

100 

147 

.626 

.674 

.791 

5570 

8576 

1  4146 

3600 

3604 

5963 

1  57 

16 

34 

200 

293 

646 

696 

761 

5670 

8576 

1  4146 

MM 

4164 

6656 

1  76 

Croit  ChwU 

34 

25 

37 

710 

543 

1  242 

6667 

1  2776 

2  1643 

1926 

3766 

5796 

301 

•‘/I.  •  .»♦« 

34 

100 

14T 

.707 

540 

1.270 

6667 

1  2776 

2  1643 

1926 

3610 

5712 

296 
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1 
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557 

1  205 
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47 
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775 
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1  110 
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1  2512 
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4292 

7303 
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TABLE  II  SUMMARY  OF  PARACHUTE  SHAPE  TEST  RESULTS 
FOR  24-OORE  AND  SO-GORE  CONFIGURATIONS 
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X.  List  of  Symbols 

Steady-state  projected  area  of  the 
Inflated  parachute,  ft^ 

■  Instantaneous  canopy  niouth  area,  ft 


-  Steady-state  Inflated  r.outh  area, 
ft^ 

-  Acceleration,  ft/sec^ 

-  Maximum  inflated  parachute  diameter 
of  ijore  malnseam,  ft 

-  Minor  axis  of  the  ellipse  bounded 
by  the  major  axis  (2a)  and  the  vent 
of  the  canopy,  ft 

-  Minor  axis  of  the  ellipse  which 
includes  the  skirt  hea  of  the 
canopy,  ft 

-  Effective  porosity 

-  Parachute  coefficient  of  drag 

-  Parachute  pressure  coefficient, 
relates  internal  and  external 
pressure  (ap)  on  canopy  surface  to 
the  dynamic  pressure  of  the  free 
stream 

-  Komlnal  diameter  of  the  aerodynamic 

decelerator  -  ft 

-  Instantaneous  force,  lbs 

-  Steady-state  drag  force  that  would 
be  produced  by  a  fully  cpen  para¬ 
chute  at  velocity  Vs,  lbs 

-  Constructed  strength  of  the  para¬ 
chute,  lbs 

-  Maximum  opening-shock  force,  lbs 

p 

-  Gravitational  acceleration,  ft/sec 

-  Permeability  constant  of  canopy 
cloth 

-  Mass,  slugs 

-  Mass  ratio  -  ratio  of  the  mass  of 
the  retarded  hardware  (including 
parachute)  to  a  mass  or  atmosphere 
contained  in  a  right  circular 
cylinder  of  length  (Vstg),  face  area 
(CjjS^),  and  density  (o) 

-  Mass  flow  ratio  -  ratio  of  atmo¬ 
sphere  flowing  through  a  unit  cloth 
area  to  the  atmosphere  flowing 
through  a  unit  inlet  area  at 
arbitrary  pressure 

-  Permeability  constant  of  canopy 
cloth 


q 

s 


S  «=A 
o  so 


V 


lo 


p 

n 


max 


S.F. 


-  Cloth  permeability  -  rate  of  air¬ 
flow  through  a  cloth  at  an  arbi¬ 
trary  differential  pressure, 

ft  Vft^/sec 

-  Dynamic  pressure,  Ic/ft' 

-  Instantaneous  inflated  canopy 
surface  area,  ft^ 

Canopy  surface  area,  ft" 

-  Instantaneous  tine,  sec 

-  Reference  time  v/hen  the  parachute 
has  reached  the  design  drag  area 
for  the  first  time,  sec 

-  Canopy  inflation  tine  when  the 
inflated  canopy  has  reached  its 
r.aximura  physical  size,  sec 

-  Air  velocity  through  cloth  In 
effective  porosity,  ft/sec 

-  Fictitious  theoretical  velocity 
used  in  effective  porosity,  ft/sec 

-  Instantaneous  systes  velocity, 

-  System  velocity  at  the  time  t  =  t  , 
ft/sec 

-  System  velocity  at  the  end  of 
suspension  line  stretch,  ft/sec 

-  Volume  of  air  which  must  be  collec¬ 
ted  during  the  inflation  process, 
ft3 

-  Hardware  weight,  lb 

-  Instantaneous  shock  factor 


-  Shock  factor  at  the  tine  t  <»  t^ 

-  Air  density,  slugs/ft^ 

-  Patio  of  parachute  projected  mouth 
area  at  line  stretch  to  the  steady- 
state  projected  area 

-  Instantaneous  elongation 

-  Maximum  elongation 

-  Initial  eloneation  at  the  beginning 
of  the  elastic  phase  of  inflation 

-  Parachute  safety  factor  =  ^c^^nax 
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Appendix  B 


A  GUIDE  FOR  THE  USE  OF  APPENDIX  A 


At  first  reading.  Appendix  A  may  appear  to  be  a  complicated  s 
tern  or  analysis  because  of  the  many  formulae  presented.  Actually 
once  understood,  the  technique  is  straightforward  and  uncorapl icat 
The  author  has  attempted  to  simplify  the  algebra  wherever  possibl 
This  appendix  presents,  in  semi -outline  form,  a  guide  to  the 
sequence  of  calculations  because  the  analysis  does  require  use  of 
formulae  from  the  text,  not  necessarily  in  the  order  in  which  the 
were  presented.  Also,  the  user  can  be  referred  to  graphs  of  per¬ 
formance  to  illustrate  effects. 


4.  W,  lb,  system  weight  (including  weight  of  the  parachute) 
from  design  requirements. 


5.  ft-^,  this  volume  of  air,  which  is  to  be  collected 
during  inflation,  is  calculated  from  the  steady-state  inflated 
shape  geometry  of  the  particular  parachute  type.  The  nomenclature 
is  described  in  Figure  24,p.A-14.  When  Dq  or  Dp  is  known,  a  can  be 
calculated  from  data  in  Table  I  and  Table  II,  p.  A- 15,  for  various 
parachute  types  and  number  of  gores.  Then  the  geometric  volume  Vo 
can  be  calculated  by  Equation  (31),  p.  a-14,  with  appropriate 
values  of  b/a  and  b'/5  from  the  tables. 

6.  Aj^q,  ft^,  steady-state  canopy  mouth  area 
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where  N/a,  b/a,  and  b'/S  are  available  from  Tables  I  and  II  for  the 
particular  type  of  parachute  and  number  of  gores. 

7.  Ago/  ft^,  canopy  surface  area  = 

8.  C^,  pressure  coefficient,  see  Figure  18,  p.  A-12.  A 
constant  Cp  =  1.7  for  all  altitudes  seems  to  yield  acceptable 
results. 


9.  Constants  k  and  n  are  derived  from  measurements  of  the  air 
flow  through  the  cloth.  Only  k  is  needed  for  Equation  (14) ,  but  n 
is  also  required  for  Equation  (13).  These  parameters  can  be  deter¬ 
mined  for  any  cloth  using  the  technique  described  beginning  on 
p.  A-12.  The  two-point  method  is  adequate  if  the  AP  across  the 
cloth  is  in  the  range  of  AP  for  actual  operation.  Check-points  of 
cloth  permeability  can  be  measured  and  compared  to  calculated 
values  to  verify  agreement.  If  the  data  are  to  be  extrapolated  to 
operational  AP's  greater  than  measured,  a  better  method  of  deter¬ 
mining  k  and  n  from  the  test  data  would  be  a  least  squares  fit 
through  many  data  points.  This  way  errors  due  to  reading  either 
of  the  two  points  are  minimized. 

11.  Step  1 

Calculate  the  reference  time  to  by  use  of  Equations  (13)  or 
(14),  p.  A-7.  If  the  deployment  altitude  is  50,000  feet  or  higher. 
Equation  0.4)  is  preferred  due  to  its  simplicity.  For  altitudes 
from  sea  level  to  50,000  feet,  Equation  (13)  is  preferred.  Figure 

12,  p.  a-8/  shows  the  effect  of  altitude  on  t^,  and  can  be  taken  as 

a  guide  for  the  user  to  decide  whether  to  use  Equation  (13)  or  (14) . 
One  should  keep  in  mind  that  the  opening  shock  force  can  be  a 
.strong  function  of  inflation  time,  so  be  as  realistic  as  possible. 

If  Equation  (13)  is  elected,  the  method  in  use  at  the  NSWC/WO  is 
to  program  Equation  (13)  to  compute  the  parachute  volume,  Vq,  for 
an  assumed  value  of  tQ.  Equation  (14) ,  because  of  its  simplicity, 
can  be  used  for  a  first  e.,timate  of  to  at  all  altitudes.  The  com¬ 
puted  canopy  volume  is  then  compared  to  the  canopy  volume  calculated 
front  the  geometry  of  the  parachute  as  per  Equation  (31),  p.  A-14. 

If  the  volvimc  computed  from  the  mass  flow  is  within  the  volume 
computed  from  the  geometry  within  plus  or  minus  a  specified  delta 
volume,  the  time  to  is  printed  out.  If  not  within  the  specified 
limits,  to  is  adjusted,  and  a  new  volume  calculated.  For  a  35-foot 
Dq,  T-10  type  canopy,  I  use  plus  or  minus  10  cubic  feet  in  the 
volume  comparison.  The  limit  would  be  reduced  for  a  parachute  of 


smaller 

If 

Vq  calculated 

0 

> 

II 

geometry  ± 

10,  then 

print  answer. 

If 

Vo  calculated 

^  '^o 

geometry  + 

10,  then 

correct  to  as  follows 

Otomatry 

*0  calculatad  (B-2) 
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The  new  value  of  is  substituted  in  the  "do  loop"  and  the 
volume  recomputed.  This  calculation  continues  until  the  required 
volume  is  within  the  specified  limits. 

III.  Calculate  to  corrected  for  initial  area.  The  t^  of  Section 
II  assumes  that  the  parachute  inflated  from  a  zero  initial  area. 

If  this  is  a  reasonable  assumption  for  the  particular  system  under 
study,  then  the  mass  ratio  can  be  determined  from  Equation  (6) , 

p.  A-'l.  For  n  =  0  if  the  value  of  M  i  0.19,  then  a  finite  state 
of  deployment  exists,  and  the  time  ratio  of  occurrence  and  the 
maximum  shock  factor  can  be  determined  from  Equations  (9)  and  (10) , 
respectively,  on  p.  A-5.  If  n  0,  then  the  limiting  mass  ratio 
for  finite  operations  will  rise  slightly  as  described  in  Appendix 
C.  Figures  C-1  and  C-2  illustrate  the  effects  of  initial  area  on 
limiting  mass  ratios  and  shock  factors  respectively.  If  the  mass 
ratio  is  greater  than  the  limiting  mass  ratio  (Ml),  then  the 
maximum  shock  force  occurs  at  a  time  greater  than  to  and  the 
elasticity  of  the  materials  must  be  considered  (see  Section  VI). 

If  n  5^  0,  then  the  reference  time,  t^,  will  be  reduced,  and 
the  mass  ratio  will  rise  due  to  partial  inflation  at  the  line 
stretch.  Figures  9  and  10,  p.  A-6,  illustrate  the  effects  of 
initial  area  on  the  velocities  and  shock  factor  during  the 
"unfolding"  inflation.  Equation  (15),  p.A-9,  can  be  used  to 
correct  to  calculated  for  the  cases  where  n  =  Ai/Ac.  If  the  initial 
value  of  drag  area  is  known,  Equation  (16) ,  p.  A-9,  can  be  used  to 
correct  tQ  and  rechecked  for  limiting  mass  ratios  versus  n  in 
Appendix  C. 

IV.  Opening  shock  calculations  in  the  elastic  phase  of  inflation. 
It  has  been  considered  that  from  time  t  =  0  to  t  =  t©  the  para¬ 
chute  has  been  inelastic.  At  the  time  t  =  to  the  applied  aero¬ 
dynamic  load  causes  the  materials  to  stretch  and  the  parachute 
canopy  increases  in  size.  The  increased  size  results  in  an 
increase  in  load,  which  causes  further  growth,  etc.  This  sequence 
of  events  continues  until  the  applied  forces  have  been  balanced  by 
the  strength  or  materials.  The  desigm-r  must  insure  that  the 
constructed  strength  of  the  materials  is  sufficient  to  resist  the 
applied  loads  for  the  material  elongation  expected.  Use  of 
materials  of  low  elongation  should  result  in  lower  opening  shock 
forces  as  Cpf^max  is  reduced . 

When  the  mass  ratio  of  the  system  is  greater  than  the  limiting 
mass  ratio,  the  elasticity  of  the  materials  and  material  strength 
determine  the  maximum  opening  shock  force.  Tl'.e  maximum  elongation 
tmax  ultimate  strength  of  the  materials  are  known  from 

tests  or  specifications.  The  technique  begins  on  p.  A-9. 

At  the  time  t  =  t^,  calculate  the  following  quantities  for  the 
particular  values  of  M  and  t,. 

a.  V^/Vg  from  Equation  (18),  p.  A- 9. 
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b.  Xq  from  Equation  (22) ,  p.  A-10. 

c.  Cq  from  Equation  (23),  p.  A-10. 

d.  Determine  CpS^jay/CoSQ  from  Figure  15,  p.  A-10. 

e.  Calculate  the  inflation  time  ratio  tf/tQ  from  Equation 
(24) ,  p.  A-10. 

f.  Calculate  the' maximum  shock  factor  from  Equation  (19), 
p .  A- 9 . 

g.  Calculate  the  opening  shock  force  where 

h.  Calculate  filling  time,  tf(sec) 

V.  In  order  to  simplify  the  required  effort,  the  work  sheets  of  Table 
B-1  are  included  on  pages  B-5  through  B-9  to  aid  the  engineer  in 
systematizing  the  analysis.  The  work  sheets  should  be  reproduced  to 
provide  additional  copies. 
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Table  B-1.  Opeiiing  Shack  Force 


CALCULATION  WORK  SHEETS 


1. 

2. 


Parachute  type  - 
System  parameters 

a.  System  weight,  W  (lb) 

2 

b.  Gravity,  g  (ft/sec  ) 

c.  Deployment  altitude  (ft) 

d.  Deployment  air  density,  p  (slugs/ft^) 

e.  Velocity  at  line  stretch,  Vg  (fps) 

f.  Steady  state  canopy  data 

(1)  Diameter,  Dq  (ft) 

(2)  Inflated  diameter,  2i  (ft); 

,  ff  p 

(3)  Surface  area,  S^  (ft'^); 

(4)  Drag  area,  (ft2);  Cp  x 

(5)  Mouth  area,  A^q  (ft^) 


(6)  Volume,  (ft^) 


^MO 


Vo 


using  technique  beginning  on  k 

Note:  Permeability,  is  usually  n 

measured  as  f t^/f t'^/min. 

For  these  calcularions 
permeability  must  be 
expressed  as  ft^/rtvaec 


*  Data  for  these  calculations  are  listed  in  Tables  1 
and  2 ,  p .  A- 1 5 . 


g.  Cloth  data 

(1)  k  1  Calculate 

(2)  n  (  p.  A-12. 


ft* 


ft* 
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Table  B-1 •  Opening  Shock  Force 

(cont 'd) 

(3)  determine  maximum  elongations  from 

pull  test  data  of  Joints,  seams,  lines,  etc.  Use 
minimum  (determined  from  tests. 

(4)  Cp;  pressure  coefficient 
h.  Steady  state  drag,  Fg  (lb). 


i.  Parachute  constructed  strength,  F^,  (lb);  deter¬ 
mined  from  data  on  efficiency  of  seams,  joints,  lines. 
Constructed  strength  is  the  minimum  load  required  to 
fall  a  member  times  the  number  of  members. 

3.  Force  calculations 

a.  Calculate  t^  for  n  ■  0;  eq.  14,  p.  A-7. 


p  PqYo  r  _ 

^ _  2W  L  A,go  -  A,o  k/CpPy  ^ 

‘'sVo  L'  ^  ^ 


Check  Figure  13,  p.  A-8,  for  advisability  of 
using  eq.  13,  p.  A-7. 

b.  If  n  ■  0,  proceed  with  steps  c  through  e. 

If  n  0,  go  to  step  f. 

c.  Mass  ratio,  M;  eq.  6,  p.  A-4 


p»V.Cd»- 


d.  If  M  <  4/21  for  n  *  0,  then  finite  mass 
deployment  is  Indicated. 

(1)  Time  ratio  at  eq.  9,  p.  A-5 


1 

j.  /iimV 


(2)  Max  shock  factor,  Xj ;  eq. ,  10,  p.  A-5 
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Table  B-1 .  Opening  Shock  Force 

(Cont'd) 

(3)  Max  shock  force, 


•  m«x.  " 


e.  If  M  >  4/21;  then  intermediate  mass  or  infinite 
mass  deployment  is  indicated  and  the  elasticity  of 
materials  is  Involved.  Calculate  Che  trajectory  con¬ 
ditions  at  time  t  *  t„. 


(1)  Velocity  ratio  @  t 


T-q  for  n 


1  + 


7M 


(2)  Shock  factor  @  t  »  for  n  -  0 

['•m]  VV 

(3)  Initial  elongation,  Cq',  eq.  23,  p.  A-10 


Cr.S_ 


A-10 


(4)  Determine  from  Figure  15,  p.  A-IO 

Vo 

(5)  Calculate  inflation  time  ratio,  tf  ;  eq.  24, 


CoSn...\* 


(6)  Calculate  maximum  shock  factor,  x^  max'* 
19,  p.  A-9  /tA* 


(cT 


(7)  Calculate  maximum  shock  force,  (lb). 
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Table  B-1 .  Opening  Shock  Force 

(cont ’d) 


(8)  Inflation  time,  sec  «  V  “  *o 


(i) 


f.  If  n  0,  correct  t^  for  initial  area  effects; 

eq.  16,  p.  A-9  ^  r  /CoSjV'^  _ 

*0  M*  )  I  *« 

g.  Mass  Ratio,  M,  eq.  6,  p.  A-4  ^  e/  J 


M 


2vy 


o 


h.  Calculate  limiting  mass  ratio, 


If  M  <;  ,  fin. .e  mass  deployment  is  indicated  and 

can  be  determined  by  eq.  8,  p.  A-5  by  assuming  values 
of  c/Cq  and  plotting  the  data  using  the  methods  of 
Appendix  C. 

i.  If  M  >  M^,  then  intermediate  mass  or 
infinite  mass  deployment  is  indicated  and  the 
elasticity  of  materials  is  involved.  Calculate 
the  trajectory  conditions  at  time  t  ■  tg. 


p.  A-9 


(1)  Velocity  ratio  @  t  ■  t-  for  n  0 ;  eq.  18, 


V. 


»  + 1  +  alLjzI  ♦ 

M  7  2 


22,  p.  A-10 


(2)  Shock  factor  5  t  -  t^,  for  n  0;  eq. 


(3)  Initial  elongation,  eq,  23,  p.  A-10 


F- 


CnS 

(4)  Determine  from  Figure  15,  p.  A-IO 

Vo 
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(5) 

p.  A- 10 


Table  B-1.  Opening  Shock  Force 
Calculate  inflation  time  ratio.  _  ;  eq-.  24, 


(6)  Calculate  maximum  shock  factor,  x. 
19,  P.  A-9  -y 

Iw 


^i«n«  r  V, 


(7)  Calculate  maximum  shock  force.  (lb) 

msK 

(8)  Calculate  inflation  time,  tffsec) 

♦f  -  *o 


(^) 
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Appendix  C 

EFFECT  OF  INITIAL  AREA  RATIO  OH  THE  LIMITING  MASS  RATIO 
AND  SHOCK  FACTOR  FOR  THE  FINITE  STATE  OF 
SOLID  CLOTH  PARACHUTE  DEPLOYMENT 


Very  low  mass  ratios  are  indicative  of  finite  mass  parachute 
deployment,  wherein  the  maximum  shock  force  occurs  before  the 
parachute  is  fully  inflated  during  the  unfolding  phase  of  deploy¬ 
ment.  As  the  mass  ratio  is  increased,  the  maximum  shock  force 
occurs  later  in  the  inflation  process.  At  some  value  of  mass 
ratio,  the  maximum  shock  force  will  occur  at  the  time  t©.  This 
particular  mass  ratio  is  defined  as  the  limiting  mass  ratio  (Ml) 
for  finite  mass  deployment.  A  further  increase  in  mass  ratio 
will  result  in  the  maximum  shock  force  occurring  after  the  para¬ 
chute  has  achieved  the  design  drag  area  (Cj^Sq)  for  the  first  time. 

Equation  (8),  from  p.  A-5,  Appendix  A,  defines  the  instantaneous 
shock  factor  during  the  unfolding  phase  of  parachute  deployment. 


This  expression  is  to  be  analyzed  for  the  following  purposes: 

a.  Determine  the  effect  of  the  initial  area  ratio  (n)  on  the 
limiting  mass  ratio  ,.Ml)  • 

b.  Determine  the  variations  of  the  instantaneous  shock  factor 
during  the  unfolding  phase  of  deployment  as  a  function  of  limiting 
mass  ratio  and  n. 

c.  Determine  the  expression  for  the  time  of  occurrence  of  the 
maximum  shock  force  for  finite  mass  ratios  less  than  the  limiting 
mass  ratio. 

The  maximum  shock  force  occurs  at  the  point  in  finite  mass 
deployment  where  dxj^/dt  =  0.  Setting  the  derivative  dx^/dt  =  0 
and  solving  for  mass  ratio  as  a  function  of  n  and  t/t^  results  in 
the  following  equality: 


C-1 
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M-^ - ^ - — _ \«/  J  j1 


t  ■ 

«•. 


(CD 


^  Since  the  limiting  mass  ratio  occurs  at  t/t^  =  1,  Equation  (C-1) 

I  can  be  reduced  to: 


M 


L 


1 

3(1  -fj) 


(C-2) 


The  effects  of  initial  area  ratio  on  the  limiting  mass  ratio  are 
described  in  Equation  (C-2)  and  Figure  C-1.  Note  that  the  time  of 
occurrence  of  the  maximum  shock  force  for  n  =  0  in  Equation  (C-1) 
is : 


'a«x, 


mM 


(C-3) 


which  is  the  same  as  Equation  (9),  p.  A-5,  Appendix  A. 

The  variation  of  the  instantaneous  shock  factor  during  the 
unfolding  phase  of  deployment  for  limiting  mass  ratios  is  presented 
in  Figure  C-2.  Initial  area  at  the  beginning  of  inflation  causes 
the  initial  force  to  increase,  but  this  is  compensated  for  by 
reduced  maximum  shock  forces.  As  n  increases,  the  initial  loads 
can  be  greater  than  the  maximum  shock  force.  However,  values  of 
n  are  usually  snail  and  depend  on  the  deployment  systems  for 
magnitude  and  repeatability.  Values  of  n  =  0.4  are  more  repre¬ 
sentative  of  a  parachute  being  disreefed  rather  than  initially 
deployed.  This  does  demonstrate,  however,  that  the  analysis  pre¬ 
sented  in  Appendix  A  can  be  adapted  to  the  disreefing  of  solid  cloth 
parachutes  by  considering  the  next  stage  to  be  a  deployment  with  a 
large  value  of  n.  Variation  in  initial  area  is  one  of  the  causes 
of  variation  in  opening  shock  forces.  The  variation  of  opening 
shock  forces  for  finite  mass,  intermediate  mas.s,  and  infinite  mass 
states  of  deployment  can  be  evaluated  by  successive  calculations 
with  various  expected  values  of  n. 

For  known  mass  ratios  less  than  the  limiting  mass  ratio,  the 
time  of  occurrence  of  the  maximum  shock  force  can  be  ascertained 
from  Equation  C-1.  If  n  approaches  zero,  then  the  time  ratio  of 
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FIGURE  C-1.  EFFECT  OF  INITIAL  AREA  RATIO  ON  THE  UMITINQ 
MASS  RATIO  FOR  THE  FINITE  STATE  OF  PARACHUTE 
DEPLOYMENT  FOR  SOLID  CLOTH  PARACHUTES 
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occurrence  of  the  maximum  shock  force  can  be  initially  estimated 
from  Equation  (C-3) ,  or  determined  by  plotting 


as  in  Figure  C-2. 
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